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Abstract 

The production of jets in low Q 2 ep scattering (photoproduction) and in low Q 2 e + e~ scatter- 
ing (77 scattering) allows for testing perturbative QCD and for measuring the proton and photon 
structure functions. This requires exact theoretical predictions for one- and two-jet cross sections. 
We describe the theoretical formalism, giving sufficient details, for calculating the direct and re- 
solved processes in jp and 77 reactions in next-to-leading order QCD. We present the complete 
analytical results for the Born terms, the virtual, and the real corrections. To separate singular and 
regular regions of phase space we use the phase space slicing method with an invariant mass cut-off. 
In this way, all soft and collinear singularities are either canceled or absorbed into the structure 
functions. Using a flexible Monte Carlo program, we evaluate the cross sections numerically and 
perform various tests and comparisons with other calculations. We consider the scale dependence 
of our results and compare them to data from the experiments HI and ZEUS at HERA and OPAL 
at LEP. 
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1 Introduction 

The analysis of jet production in various high energy processes has become a major field for testing 
perturbative QCD. Recently jet production in 77 processes, where the two photons of very small 
virtuality are produced in the collision of electrons and positrons, has come into focus after data have 
been collected at the TRISTAN Q and LEP [0j colliders. In addition jet production in deep inelastic 
(high Q 2 ) |J and low Q 2 ep scattering (equivalent to photoproduction) has been measured at the two 
HERA experiments HI f§ and ZEUS §. 

Jet production in 77 and 7p collisions has several similarities. At very small Q 2 , where q 
(Q 2 = —q 2 ) is the four-momentum transfer of the electron (positron) producing the virtual photons 
in the 77 or jp initial state, the emission of the photon can be described in the Equivalent Photon 
Approximation. The spectrum of the virtual photons is approximated by the Weizsacker- Williams 
(WWA) formula, which depends only on y = E^/E e , the fraction of the initial electron (positron) 
energy E e transferred to the photon with energy E^, and on Qj^ax ( or $max)> which is the maximal 
virtuality (or the maximal electron scattering angle) allowed in the experimental set-up. 

Concerning the hard scattering reactions both processes have so-called direct and resolved com- 
ponents. Thus, in leading order (LO) QCD the cross section (1(77 — ► jets) receives contributions 
from three distinct parts: (i) the direct contribution (DD), in which the two photons couple directly 
to quarks, (ii) the single-resolved contribution (DR), where one of the photons interacts with the 
partonic constituents of the other photon, (iii) the double-resolved contribution (RR), where both 
photons are resolved into partonic constituents before the hard scattering subprocess takes place. In 
the DD component (in LO) we have only two high-py jets in the final state and no additional spectator 
jets. In the DR contribution one spectator jet originating from low transverse momentum fragments 
of one of the photons is present, and in the RR component we have two such spectator or photon 
remnant jets. In the case of 7p collisions one of the photons is replaced by the proton which has no 
direct interaction. Then we have only the DR and RR components which are usually referred to as 
the direct and the resolved contribution. This means that the DR component for 77 — » jets has the 
same structure as the direct contribution of 7^ — > jets and the RR part for 77 — > jets is calculated in 
the same way as the resolved cross section for 'yp — > jets. 

Of course, to calculate the resolved cross sections we need a description of the partonic constituents 
of the photon. These parton distributions of the photon are partly perturbative and non-perturbative 
quantities. To fix the non-perturbative part one needs information from experimental measurements 
to determine the fractional momentum dependence at a reference scale Mq . The change with the fac- 
torization scale M 2 < Mq is obtained from perturbative evolution equations. Most of the information 
on the parton distribution functions (PDF) of the photon comes from photon structure function (i 7 ^) 
measurements in deep inelastic ej scattering, where, however, mainly the quark distribution function 
can be determined. In 77 and 7p high-py jet production also the gluon distribution of the photon 
enters which can be constrained by these processes. 

When we want to proceed to next-to-leading order (NLO) QCD the following steps must be taken: 
(i) The hard scattering cross section for the direct and resolved photon processes are calculated up to 
NLO. (ii) NLO constructions for the PDF's of the proton and the photon are used and are evolved 
in NLO up to the chosen factorization scale via the Altarelli-Parisi equations and convoluted with 
the NLO hard-scattering cross sections, (iii) To calculate jet cross sections we must choose a jet 
definition which may be either a cluster or a cone algorithm in accordance with the choice made in 
the experimental analysis. 

There exist several methods for calculating NLO corrections of jet cross sections in high energy 
reactions §. As in our previous work we apply the phase space slicing method with invariant mass 
slicing to separate infrared and collinear singular phase space regions. In this approach the contri- 
butions from the singular regions are calculated analytically with the approximation that terms of 
the order of the slicing cut are neglected. In the non-singular phase space regions the cross section is 



obtained numerically. This method allows the application of different clustering procedures, the use 
of different variables for describing the final state together with cuts on these variables as given by the 
measurement of the jet cross sections. This method has been used for the calculation of the NLO DD 
and DR cross sections in the case of 77 — ► jets J7[ and for the calculation of the direct cross section 
in the case of jp — ► jets M. It is obvious that the NLO calculation of the DR cross section for the 77 
case is the same as of the NLO direct cross section for the 7p case. 

In the calculation of the 77 DD and DR or 7^ direct cross section one encounters the photon- 
quark collinear singularity. This is subtracted and absorbed into the quark distribution of the photon 
in accord with the factorization theorem. This subtraction at the factorization scale M 2 produces an 
interdependence of the three components in the 77 and the two components in the 7p reaction, so that 
a unique separation in DD, DR, and RR (for 77) and in direct and resolved (for 'jp) contributions is 
valid only in LO. The calculation of these subtraction terms and of the contributions from the other 
singular regions has been presented with sufficient details in ffl, M. The phase space slicing method was 
applied also for the calculation of the NLO correction of the 77 RR contribution Q and the resolved 
7P cross section [1C]. Results for the inclusive two-jet cross sections together with a comparison to 



recent experimental data were presented for 77 — ► jets in [11| and for 7J) — » jets in [jT2| , respectively. 
In these two papers the specific calculations of the corresponding resolved cross sections in NLO 
were not explicitly outlined, in particular it was not shown, how the different infrared and collinear 
singularities cancel between virtual and real contributions and after the subtraction of collinear initial 
state singularities. In this work we want to fill this gap. We describe the analytic calculation of the 
various resolved terms in the specific singular regions needed for the NLO corrections. Furthermore 
we check the analytical results by comparing with results obtained with other methods. In order to 
have a complete presentation we include also the calculation of the NLO correction to the 77 DD 
cross section already presented in j/J. We also include the details of the DR contributions in the form 
of the direct jp cross section taken from ||]. Details of other material, as for example, various jet 
definitions, which were only mentioned in our previous work and which now become more and more 
relevant in the analysis of the experimental data, will also be presented. Since the calculation of the 
RR (DR) cross section for 77 — » jets is the same as the calculation of the jp — > jets resolved (direct) 
cross section we concentrate on the jp case when we present the details of the calculations. So for the 
77 case we present only the DD contribution which has no analogy in the 7^ case. We come back to 
the 77 case when we show the numerical results for specific cases including all three components. 

We organize this work in seven main sections: After this introduction we relate the experimental 
ep scattering to photon parton scattering in section 2. We describe the Weizsacker- Williams approxi- 
mation, discuss the proton and photon PDF's and explain the experimental and theoretical properties 
of various jet definitions. Furthermore, section 2 contains the master formulas for one- and two-jet 
cross sections. These will be calculated in section 3 in LO and in section 4 in NLO. In both sections, 
we calculate the relevant phase space for 2^2 and 2^3 scattering, respectively. The Born matrix 
elements are contained in section 3. In section 4 we present the virtual one-loop matrix elements 
and the tree-level 2 — ► 3 matrix elements, which are then integrated over singular regions of phase 
space. Next, we demonstrate how all ultraviolet and infrared poles in the NLO calculation cancel or 
are removed through renormalization and factorization into PDF's. A detailed numerical evaluation 
of jet cross sections in 7p scattering with the purpose to compare with results of other work and to 
make consistency checks is contained in section 5. We study the renormalization and factorization 
scale dependence of one- and two-jet cross sections including the direct and resolved components. We 
also show results for some specific inclusive one- and two-jet cross sections and compare them with 
experimental data, in case they are available, to demonstrate the usefulness of our methods. Section 
6 contains the corresponding numerical studies and comparisons to data for 77 scattering. The final 
conclusions are left for section 7. 



2 Photoproduction of Jets at HERA 

In this section we set up the general framework for theoretical predictions of the photoproduction of 
jets in electron-proton scattering. This includes the separation of the perturbatively calculable parts 
from the non-perturbative parts of the cross section and linking the electron-proton to photon-proton 
scattering. This link will be discussed first in section 2.1 and consists in the Weizsacker- Williams or 
Equivalent Photon Approximation. 

The framework of the QCD improved parton model for protons will shortly be discussed in section 
2.2. This is necessary, since protons are not pointlike but composed of three valence quarks and sea 
quarks and gluons. Perturbative QCD is not applicable at distances comparable to the size of the 
proton, but only at small partonic scales due to the asymptotic freedom of QCD. Therefore, the parton 
content of the proton is the domain of non-perturbative QCD and has to be described with universal 
distribution functions. The scale dependence of these functions is governed by the Altarelli-Parisi 
equations. 

A similar concept applies for resolved photons, which are discussed in section 2.3. Contrary to di- 
rect photons, which are obviously pointlike, resolved photons can be considered to have a complicated 
hadronic structure like protons. However, they have a different valence structure than protons. Fur- 
thermore, a complex relationship between direct and resolved photons will show up in next-to-leading 
order of QCD. 

Having related the initial state particles electron and proton in the experiment to the photons 
and partons in perturbative QCD, we can turn our attention to the final state particles. Section 2.4 
describes how the interpretation of partons as jets changes from leading to next-to-leading order of 
QCD. Several jet definition schemes are discussed with respect to their theoretical and experimental 
behavior. 

The last section 2.5 summarizes the different ingredients of the calculation and contains the master 
formula? for one- and two-jet photoproduction. Also, the numerous analytical contributions in leading 
and next-to-leading order of QCD are organized in a tabular form for transparency. 

2.1 Photon Spectrum in the Electron 

In photoproduction, one would like to study the hard scattering of real photon beams off nuclear 
targets. This has been done in fixed target experiments, like NA14 at CERN [l^] or E683 at Fermilab 
|fL4| l, where real photons are produced e.g. in pion decay with energies of up to 400 GeV j[j|]. If higher 
energies are required, one must resort to spacelike, almost real photons radiated from electrons. This 
method is employed at the electron-proton collider HERA. There, electrons of energy E e = 26.7 GeV 
and recently positrons of energy E e = 27.5 GeV produce photons with small virtuality Q 2 , which then 
collide with a proton beam of energy E p = 820 GeV. This corresponds to photon energies of up to 50 
TeV in fixed target experiments. 

On the theoretical side, the calculation of the electron-proton cross section can be considerably 
simplified by using the Weizsacker- Williams or Equivalent Photon Approximation. Here, one uses 
current conservation and the small photon virtuality to factorize the electron-proton cross section into 
a broad-band photon spectrum in the electron and the hard photon-proton cross section. Already 
in 1924, Fermi discovered the equivalence between the perturbation of distant atoms by the field 



of charged particles flying by and the one due to incident electromagnetic radiation [16]. His semi 



classical treatment was then extended to high-energy electrodynamics in 1933-1935 by Weizsacker 



[ p"7| ] and Williams [18 1 independently, who used a Fourier analysis to unravel the predominance of 
transverse over longitudinal photons radiated from a relativistic charged particle. In the fifties, Curtis 
[ |l9| ] and Dalitz and Yennie |20[] gave the first field-theoretical derivations and applied the approximation 
to meson production in electron- nucleon collisions. Chen and Zerwas used infinite-momentum- frame 
techniques for an extension to photon bremsstrahlung and photon splitting processes [21]. For a recent 



review of the various approximations see [22] and for the application to 77 processes see [23]. 



Let us consider the electro-production process 

electron(fc) + proton(p) — ► electron(fc') + X (2-1) 

as shown in figure ||, where k, k', and p are the four-momenta of the incoming and outgoing electron 



P 




>X 



Figure 1: Electron-proton scattering with single photon exchange. 
and the proton, respectively. X denotes a generic hadronic system not specified here, q = k — k' is the 



momentum transfer of the electron to the photon with virtuality Q^ 



0, and the center-of-mass 



energy of the process is ^[s~h = \/(k +p) 2 = 295.9 GeV (= 300.3 GeV for the positron beam). We 
restrict ourselves to one-photon exchange without electroweak Z° admixture. Two-photon exchange 
is suppressed by an additional order of a, and Z° exchange is suppressed by the pole in the Z° mass 
l/m| . 



The cross section of this process is 



da ep (ep -y eX) = 



1 e 2 W^ v T llv d 3 k 
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where W^ v and T^ are the usual hadron and lepton tensors. Exploiting current conservation and the 
small photon virtuality, we find 



W^T^ = AW 1 (Q 2 = 0,q.p) 
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where 

q.p k'.p 

Xa = T~ = 1 ~ T— 
k.p k.p 

is the fraction of longitudinal momentum carried by the photon. Next, we have to calculate the phase 
space for the scattered electron. Using 



k' = {E' e , 0,E' e /3' sin 6,E' e /3' cos 6) , k = (E e ,0,0,E e f3), 

P' = \l 



™ 2 e R- h _™± 



(2.5) 

(2.6) 



and integrating over the azimuthal angle, we find 

^ = 2ir(3' 2 E' e dk'dcos8 

= Tldq 2 dx a . (2.7) 

Combining eq. (|2.3| ) and ( |2.7| ) and integrating over the photon virtuality, we can factorize the electron- 
proton cross section into 

i 



da ep (ep -> eX) = I dx a F^/ e (x a )da^ p (-fp -> X), (2. 
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is the renowned Weizsacker- Williams approximation and where 

^TP^^)— agg^-^g -*-- (2.10) 

is the photon-proton cross section. We will only use the leading logarithmic contribution and neglect 
the second term in eq. ( |2.9D , also calculated by Frixione et al. [24|. 

There exist two fundamentally different experimental situations. At HERA, the scattered electron 
is anti-tagged and must disappear into the beam pipe. For such small scattering angles of the electron, 
the integration bounds Q m - m and Q max can be calculated from the equation 

H^- 
where A = E e {l + (3)(1 — x a ). Using the minimum scattering angle 8 = 0, we obtain 



Q 2 = ^^ + Ee(l + P)(A*- mi j Q2 + o{0% 



QLn = f^S (2-12) 

l-X a 

which leads us to the final form of the Weizsacker- Williams approximation used here: 

At HI and ZEUS, the maximum virtualities of the photon are given directly as Q^ax = 0-01 GeV 
and 4 GeV , respectively. If only the maximum scattering angle of the electron is known, we obtain 

Vmax = &e\* ~ x a)"max (2.14) 

with 6> max being of the order of 5°. This form of Q^ax wm be used for the calculation of the 77 cross 
sections in accordance with experimental choices of # max a t LEP. Except in equation (2.13), we will 
assume all particles to be massless in this paper. Therefore, all results are only valid in the high-energy 
limit. 

When no information about the scattered electron is available, one has to integrate over the whole 
phase space thus allowing large transverse momenta and endangering the factorization property of 
the cross section. Then, one only knows that the invariant mass of the produced hadronic system X 
has to be bounded from below, e.g. by minimal transverse momenta or heavy quark mass thresholds, 
which constrains the maximum virtuality of the photon only weakly. Possible choices are 



^;max — o ' 9 ' ' ' \ xo J 
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At e + e~ colliders, bremsstrahlung of electrons is not the only source of almost real photons. The 
particles in one bunch experience rapid acceleration when they enter the electromagnetic field of 
the opposite bunch producing beamstrahlung which depends sensitively on the machine parameters. 
Even higher luminosities can be achieved by colliding the electron beam at some distance from the 
interaction point with a dense laser beam [f25[. 



2.2 Parton Distributions in the Proton 

The first evidence that the proton has a substructure came from deep inelastic electron-proton scatter- 
ing ep — > eX. Together with muon- and neutrino-scattering, this process provides us with information 
on the distribution of partons (quarks and gluons) in the proton. In the parton model, scattering off 
hadrons is interpreted as an incoherent superposition of scattering off massless and pointlike partons. 
This is shown in figure 0, where pb is the four-momentum of the scattered parton b. If the scale Mb, 
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Figure 2: Electron-proton scattering in the parton model. 

at which the proton is probed, is large enough, the transverse momentum of the parton will be small 
and its phase space can be described by a single variable 

x 6 = ^€[0,l], (2.16) 

qp 

the longitudinal momentum fraction of the proton carried by the parton. The other partons in the 
proton will then not notice the interaction but form a so-called remnant or spectator jet. Multiple 
scattering is suppressed by factors of 1/Mjf, and only the leading twist contributes for large Mjj. 
The proton remnant carries the four- momentum p r , which again depends only on Xb and the proton 
momentum p. The proton remnant will not be counted as a jet in the following. 
We can now factorize the photon-proton cross section into 

i 
da 7 p(7p — ► jets + remnant) = ^J / &XbF b / p (xb } M b )do~ 1 b{lb — ► jets), (2-17) 

b o 

where F b i p {xb, Mjj) is the probability of finding a parton b (where b may be a quark, an antiquark, or 
a gluon) within the proton carrying a fraction Xb of its momentum when probed at the hard scale M&. 

6 



The usual choice for deep inelastic scattering, M& = Q, is not possible for photoproduction, as Q 2 ~ 
here. Instead, one takes M\, = S,Et with Ex being the transverse energy of the outgoing parton or 
observed jet and £ being of order 1. Contrary to the hard photon-parton scattering cross section dcr 7 b, 
the parton densities F b / p (xb, M 2 ) are not calculable in perturbative QCD and have to be taken from 
experiment. 

In the leading twist approximation, the parton densities extracted from deep inelastic scattering 
can be used for any other process with incoming nucleons like photoproduction or proton-antiproton 
scattering - they are universal. As stated above, they have to be determined by experiment at some 
scale Q 2 = Qq, where one parametrizes the x dependence before evolving up to any value of Q 2 
according to the Altarelli-Parisi (AP) equations [p6[| . The input parameters are then determined by 
a global fit to the data. The u, d, and s quark densities are rather well known today from muon 
scattering (BCDMS, NMC, E665) and neutrino scattering (CCFR) experiments, and the c quark is 
constrained by open charm production at EMC. The small- x region, where the gluon and the sea 
quarks become important, could, however, not be studied by these experiments. This is mainly the 
domain of the HERA experiments HI and ZEUS in DIS as well as in photoproduction. 

There exist many different sets of parton density functions. Most of them are easily accessible in 
the CERN library PDFLIB kept up to date by Plothow-Besch p?]]. We will now briefly compare the 
three parametrizations mainly in use today: CTEQ |g8|, MRS [pSj ], and GRV |3(J]. The first two are 
very similar to each other, and both take Qq = 4 GeV . Very recently, CTEQ and MRS presented new 
parametrizations (set 4 and set R, respectively) including new HERA and TEVATRON inclusive jet 
data |3l| , p2f . The GRV approach is rather different from the two mentioned above. In order to avoid 
any free additional parameters, their original version was based on the assumption that all gluon and 
sea distributions are generated dynamically from measured valence quark densities at a very low scale 
Qo ~ C(A). The QCD coupling scale A is fitted to the data and ranges from 200 to 344 MeV for the 
parametrizations considered here and for four quark flavors. 



2.3 Parton Distributions in the Photon 

Although the photon is the fundamental gauge boson of quantum electrodynamics (QED), which is 
the most accurately tested field theory, many reactions involving photons are much less well under- 
stood. This is due to the fact that the photon can fluctuate into qq pairs which in turn evolve into a 
complicated hadronic structure. At HERA, the photon radiated from the electron can thus interact 
either directly with a parton in the proton (direct component) or act as a hadronic source of partons 
which collide with the partons in the proton (resolved component, see figure |3|). In the latter case, 
one does not test the proton structure alone but also the photon structure. 




a) 



b) 



Figure 3: Generic diagrams for a) direct and b) resolved photoproduction. 



Using the factorization theorem, the cross section for photon-parton scattering can be written as 

i 
dcr 7f) (7& -> jets) = ^2 dy a F a/ ^(y a , M 2 )da ab (ab -> jets). (2.18) 

a J Q 

Here, F a /^(y a , M%) stands for the probability of finding a parton a with momentum fraction y a in the 
photon, which has to be universal in all processes. The scale M a is a measure for the hardness of the 
parton-parton cross section da a b calculable in perturbative QCD and is again taken to be of 0(Et)- 
The particle a can also be a direct photon. Then, F^i^(y a ,M 2 ) is simply given by the 5-function 
5(1 — y a ) and does not depend on the hard scale M a . 

Before HERA started taking data, information on the hadronic structure of the photon came 
almost exclusively from deep inelastic 7*7 scattering at e + e~ colliders. Similarly to deep inelastic ep 
scattering, this is a totally inclusive process well suited to define a photon structure function. Using 
y = Q 2 /(xbsh), where Q 2 denotes the virtuality of the probing photon 7*, and replacing F\ by the 
longitudinal structure function Fl(xb,Q 2 ) = F2(xb,Q 2 ) — 2xbF\(xb,Q 2 ), we can write the deep 
inelastic scattering (DIS) cross section in the following form 



d a 2ira 2 SH 



dx B dy Q 4 



{ [l + (1 - y) 2 ] F2(x B , Q 2 ) - y 2 Fl{x B , Q 2 )} , (2.19) 



where in LO F^ixBiQ 2 ) is related to the singlet quark parton density in the photon similarly as in 
deep inelastic ep scattering 

F^(x B , Q 2 ) = E x B e 2 q {F qh {x B , Q 2 ) + F- qh {x B , Q 2 )). (2.20) 

Fitting the input parameters to data is not as easy as in the proton case. First, no momentum sum 
rule applies for the photonic parton densities as they are all of LO in QED. This and the subleading 
nature of the gluonic process makes a determination of the gluon from F^ very difficult. However, a 
momentum sum rule exists for mesons so that the VMD part of the gluon is constrained. Second, the 
cross section for F^ is quite small and falls rapidly with increasing Q 2 (see eq. ( 2.1 9|) ) leading to large 
statistical errors. 

Presently, over 20 different parton distributions exist for the photon, and most of them are available 
in the PDFLIB p7|. The VMD input is insufficient to fit the data at higher Q 2 offering two alternatives: 
GRV |3*|], AFG [p4| , and SaS [35] follow the same "dynamical" philosophy as in the GRV proton case 
starting from a simple valence-like input at a low scale Qq = 0.25 ... 0.36 GeV . On the other hand, 
LAC [36] and GS [37] take a larger value for Qq = 4 ... 5.3 GeV 2 , assume a more complicated ansatz 
there, and fit the free parameters to -F^ data. LAC intended to demonstrate the poor constraints on 
the gluon assuming a very soft gluon (fits 1 and 2) as well as a very hard one (fit 3). However, LAC 
3 is ruled out now by recent e + e~ and HERA data. In their latest parametrization, GS lowered their 
input scale to Qq = 3 GeV , included all available data on i 7 ^, and constrained the gluon from jet 



production at TRISTAN |§. 



The distinction between direct and resolved photon is only meaningful in LO of perturbation 
theory. In NLO, collinear singularities arise from the photon initial state that have to be absorbed 
into the photon structure function (cf. sections 4.2.5 and 4.2.7) and produce a factorization scheme 
dependence as in the proton case. If one requires approximately the same F^[ in LO and NLO, the 
quark distributions in the MS-scheme have quite different shapes in LO and NLO. This is not the 
case if the DIS 7 -scheme of GRV is used, where the direct-photon contribution to F^ is absorbed 
into the photonic quark distributions. This allows for perturbative stability between LO and NLO 
results. Therefore, the separation between direct and resolved process is an artifact of finite order 
perturbation theory and depends on the factorization scheme and scale M a . Experimentally, one tries 
to get a handle on this by introducing kinematical cuts, e.g. on the photon energy fraction taking part 
in the hard cross section. 



2.4 Jet Definitions 

Due to the confinement of color charge, neither incoming nor outgoing partons can be observed 
directly in experiment, but rather transform into colorless hadrons. This transformation is a long- 
distance process and is not calculable in perturbative QCD. For the incoming particles, we have already 
described how the ignorance of universal parton distributions in hadrons is parametrized. A similar 
method can be used for the final state partons, employing so-called fragmentation functions for the 



inclusive production of single hadrons |3S]. The transformation of partons into individual hadrons 



forming part of the total final state was first studied by [4C]. Alternatively, one can observe beams 
of many hadrons going approximately into the same direction without the need to specify individual 
hadrons. The hadrons are then combined into so-called jets by cluster algorithms, where one starts 
from an initial cluster in phase space and ends at stable fixed points for the jet coordinates. These jet 



definitions should fulfill a number of important properties. They should [41] 

• be simple to implement in the experimental analysis, 

• be simple to implement in the theoretical calculation, 

• be defined at any order of perturbation theory, 

• yield a finite cross section at any order of perturbation theory, 

• yield a cross section that is relatively insensitive to hadronization. 

Although, in principle, hadronization of the final state should be factorizable from the hard cross 
section and the initial state, jets do indeed look quite different in e + e~ and in, at least partly, hadronic 
collisions like ep scattering. This can be attributed to the "underlying event" of remnant jet production 
from the initial state, which can interfere with the hard jets in the final state. The main problem here 
is the determination of the true jet energy and the subtraction of the remnant pedestal. 

In LO QCD, there is a one-to-one correspondence between partons and jets. This results in a 
complete insensitivity of theory to the experimentally used algorithm or to the resolution parameters. 
The experimental results depend, however, on these parameters as well as on detector properties and 
have to be corrected stepwise from detector level to hadron level to parton level. The situation can 
only be improved by going to NLO in perturbation theory. Here, the emission of one additional soft 
or collinear parton is calculable with correct treatment of the occurring singularities. A hadron jet 
can then consist not only of one, but also of two partons. It acquires a certain substructure, and will 
depend on the experimental algorithm and resolution. 

Historically, the first resolution criterion was proposed by Sterman and Weinberg |^]. It adds a 
particle to a jet if its energy is smaller than eM or its angle with the jet is less than 25, which provides 
a close link to the radiation of secondary partons. The energy cut handles the soft divergencies and 
the angular cut the collinear divergencies. There is some freedom in the choice for M. Usually, one 
takes the hard scale of the process, e.g. Q in e + e~ annihilation. 

The PETRA experiments used this (e, 5) algorithm and the so-called JADE cluster algorithm 
p3], which has the advantage of being invariant. Two particles are combined into a cluster, if their 



invariant mass Sy = (pi + Pj) is smaller than yM [44], where M is again a typical scale and y 



is of 0(1O -2 ). In this way, the soft and collinear divergencies can be described by a single cut-off 
y. There exist a number of different schemes for this algorithm regarding the invariant mass of the 
combination of the two particles and the combination of the particle four- momenta (JADE-, E-, E0-, 
P-, and PO-scheme). 

At hadron colliders, cluster algorithms tend to include hadrons from the remnant jet, which is not 
present in e + e~ collisions, into the current jet. Therefore, algorithms are preferred here which use a 
cone in rapidity-azimuth (77 — 4>) space, quite similar to the (5-condition of Sterman and Weinberg. The 
(pseudo-)rapidity -q = — ln[tan(0/2)] parametrizes the polar angle 6 between the hard jet and the beam 



axis, and 4> is the azimuthal angle of the jet around the beam axis. In the case of cone algorithms, only 
inclusive cross sections are infrared safe, where an arbitrary number of particles outside the jet cone 
can be radiated as long as they are softer than the observed jets. According to the standardization 
of the Snowmass meeting in 1990, calorimeter cells or partons i may have a distance Ri from the jet 
center provided that 

Ri = Jijh - m) 2 + (& - <Pj) 2 < R, (2.21) 



where rji and <j>i are the coordinates of the parton or the center of the calorimeter cell [^] . Typical 
values for the resolution parameter R range from 0.7 to 1, where the effects of hadronization and the 
underlying event are minimized. The transverse energy of the combined jet Et, is calculated from 
the sum of the particle E^ 

E Tj = J2 E n, (2-22) 

Ri<R 

and the jet axis is defined by the weighted averages 



1 



VJ 



Etj 
1 



Rt<R 



^TiVi 



E 



Tj 



E T A 



(2.23) 
(2.24) 



Rt<R 



In perturbative QCD, the final state consists of a limited number of partons. For a single isolated 
parton i, the partonic and jet parameters agree ((ET z ,r]i,4>i) = (Etj,i]j,(J)j)) as shown in figure |||a), 
whereas two neighboring partons i, j will form a combined jet as shown in figure [|b). Equation ( 2.21| ) 
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Figure 4: Jet cone definition according to the Snowmass convention for a) a single parton, b) two 
combined partons with distance R from the jet axis, and c) two single partons. 

only defines the distances Ri, Rj of each parton from the jet axis so that the two partons may be 
separated from each other by 



Rij = J(vi-Vj) 2 + 



- ^? < Et : +ETj , r- 

V3> - max(£ Ti ,£ T ) 



(2.25) 



If both partons have equal transverse energy, they may then be separated by as much as 2R. As 
parton j does not lie inside a cone of radius R around parton i and vice versa, one might with some 
justification also count the two partons separately as shown in figure |]c). If one wants to study only 
the highest-Er jet, this "double counting" must be excluded. 

The selection of the initiating cluster, before a cone is introduced ("seed-finding"), is not fixed by 
the Snowmass convention, and different approaches are possible. The ZEUS collaboration at HERA 
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uses two different cone algorithms: EUCELL takes the calorimeter cells in a window in rj — (f> space as 
seeds to find a cone with the highest Et- The cells in this cone are then removed, and the algorithm 
continues. On the other hand, PUCELL was adapted from CDF and starts with single calorimeter 
cells. It then iterates cones around each of them, until the set of enclosed cells is stable. In this case 
it may happen that two stable jets overlap. If the overlapping transverse energy amounts to a large 
fraction of the jets, they are merged, otherwise the overlapping energy is split. Alternatively, the 
overlap could be attributed to the nearest, to the largest, or to both jets. The question of overlapping 
jets cannot be addressed in a next-to-leading order calculation of photoproduction. There, we only 
have up to three partons in the final state, which can form at most one recombined jet and no overlap. 
Experimentally, jets of type b) in figure |] are hard to find because of the missing seed in the jet 
center. This is a problem in particular with the PUCELL algorithm, which relies on initial clusters and 
does indeed find smaller cross sections than the less affected EUCELL algorithm jl(]]. One possibility 
to model this theoretically is to introduce an additional parameter -R SC p P§] to restrict the distance of 
two partons from each other: 



Rij < min 



Et- + Et 

f 771 771 \ rt ' "Sep 

max(£r. , E Tj ) 



(2.26) 



R sep = 2R means no restriction. In figure ||, we can see that for two partons of similar or equal 




Figure 5: Jet cone definition with an additional parameter for parton-parton separation R scp : a) two 
partons with similar or equal transverse energies Et, b) two partons with large Et imbalance. 



transverse energies Et, R S ep is the limiting parameter, whereas it is the parton-jet distance R for two 
partons with large Et imbalance. Numerical studies of the R and R sep dependences will be given in 



sections 5.3 and 5.5 [46]. 

The JADE algorithm clusters soft particles regardless of how far apart in angle they are, because 
only the invariant mass of the clusters is used. This is improved in the hr or Durham clustering 



algorithm of Catani et al. [47|, where one defines the closeness of two particles by 

\2e>2 



di 



mm(E Ti ,E Tj YRf r 



(2.27) 



Rij is again defined in (n — <j>) space. For small opening angles, Rij <C 1, eq. (|2.27[) reduces to 



mm{E Tl ,E Tj ) 2 Rlj s* min(£;, Ej) 2 A6 2 ~ k T , (2.28) 



the relative transverse momentum squared of the two particles in the jet. Similarly, one can define a 
closeness to the remnant jet particles b, if they are present: 



dib — E Ti R ib . 



(2.29) 
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Particles are clustered into jets as long as they are closer than 

d {ijM < d 2 = E 2 Tt R\ (2.30) 

where R is an adjustable parameter of 0(1) analogous to the cone size parameter in the cone algorithm. 
Also, one chooses the same recombination scheme. Consequently, the fcy algorithm produces jets that 
are very similar to those produced by the cone algorithm with R = R sep . 

2.5 Hard Photoproduction Cross Sections 

We have now established the links between the experimentally observed initial and final states at 
HERA (electron, proton, and jets) and the partonic subprocess calculable in perturbative QCD. For- 
mally, we can combine eqs. (|2.8|) , ( |2.17| ), and ( |2.18 ) into 

da ep (ep — > e + jets + remnants) = (2-31) 

l l l 

= J2 J dx a F y/e (x a ) I dy a F a/l (y a , Ml) J dx b F b/p (x b , M^da^' (ab -> jets), 

a ' b 

where x a , y a , and x b denote the longitudinal momentum fractions of the photon in the electron, the 
parton in the photon, and the parton in the proton, respectively. From now on, we will use the variable 
x a as the variable for the parton in the electron with the consequence that the incoming partons have 
momenta p a = x a k and p b = x b p and eq. (|2.31|) becomes 



da ep (ep — ► e + jets + remnants) = (2.32) 

l l l 

= J2 Ux a F l/e {^) f^F ah (y a ,M 2 a ) f dx b F b/p (x b ,M 2 )da^ (ab - jets). 
*,*> o Va I Va i 

Next, one has to fix the kinematics for the photoproduction of jets. All particles are consid- 
ered to be massless. In the HERA laboratory system, the positive z-axis is taken along the pro- 
ton direction such that k = E e (l,0, 0, — 1) and p = £^(1,0,0,1). For the hard jets, we choose 
the decomposition of four-momenta into transverse energies, rapidities, and azimuthal angles pt = 
E^ (cosh rji, cos (j>i, sin (j)i,sinh.r]i). The boost from the HERA laboratory system into the ep center-of- 
mass system is then simply mediated by a shift in rapidity ?7boost = 1/2 \n(E e /E p ). Through energy and 
momentum conservation p a + p b = J2i Pi > t ne final partons are related kinematically to the invariant 
scaling variables x a and x b 

x ° = ^E^"*' (2.33) 



*b = 2E y Z E T< e ' H - ( 2 - 34 ) 



1 

The cross section for the production of an n-parton final state from two initial partons a, b, 

d^Coft-jets) = ^pWpdPSW, (2.35) 

depends on the flux factor l/(2s), where s = x a x b sn is the partonic center-of-mass energy, the n- 
particle phase space 

dPSW = J(2n) d f[^0l5 d ( Pa + Pb -±p] , (2.36) 
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where d = 4 — 2e is the space-time dimension, and on the squared matrix element |A4| 2 . The latter is 
averaged over initial and summed over final spin and color states. Since we study almost-real photons 
with Q 2 ~ and treat quarks and gluons as massless, only transverse polarizations are possible 
giving spin average factors of 1/2 for photons and quarks and 1/(2(1 — e)) for gluons in dimensional 
regularization. The quarks form the fundamental representation of the SU(3) color symmetry group 
and can therefore carry Nc = 3 different colors. The gluons belong to the adjoint representation of 
dimension 2NcCf = Nq — 1 = 8. This is the defining equation for the so-called color factor Cf- We 
therefore average the initial colors with factors of 1/Nq for quarks and 1/(2NcCf) for gluons. 

In leading order QCD, the phase space for two partons dPS^ ' in the final state and the Born 
matrix elements T have to be computed as displayed in the first line of table |l[ In next-to-leading 
order QCD, new contributions arise. Whereas the virtual corrections V also have a two-particle phase 
space (second line in table ||), the real corrections of the final state F, of the photon initial state /, 
and of the proton initial state J require the inclusion of a third outgoing parton into the phase space 
dPS^ ' (lines three to five in table [l]). We will distinguish between direct and resolved photoproduction 
in the subsequent chapters and treat the two contributions in a completely parallel way. 



Order 


Phase Space dPS^ 


Direct Matrix Element \M \ 2 


Resolved Matrix Element \M \ 2 


LO 


dPS (2) , sect. 3.1 


T 7 (,^i2, sect. 3.2 


T a b^i2, sect. 3.3 


NLO 


Vyb^i2, sect. 4.1.1 


V a b^i2, sect. 4.1.2 


dPSW, sect. 4.2.1 


^76^123, sect. 4.2.2 


F a b->123, sect. 4.2.3 


dPS (3) , sect. 4.2.4 


I-yb~>i23, sect. 4.2.5 


Iab->123, sect. 4.2.7 


^7b^i23> sect. 4.2.6 


Jab^i23, sect. 4.2.8 



Table 1: Summary of phase space parametrizations and matrix elements needed for this NLO calcula- 
tion of direct and resolved photoproduction. 

The two final state partons produced in a LO process have to balance their transverse energies 

(2.37) 
(2.38) 



Et ± = Et 2 = Et, so that relations (|2.33| ) and ( j2.34|) simplify to 

Ef 



j a 



-i n 



x h 



2E e 

Ef 
2 E n 



+ e~^) 



(e" 1 + e 



1/2 ^ 



The rapidity 772 of the second jet is kinematically fixed by Ef, i]i, and x a through the relation 



m 



In 



( 2x a E e 
\ E T 



->n 



(2.39) 



In the HERA experiments, x a is restricted to a fixed interval x ajm i n < x a < £ ajmax < 1- We shall 

disregard this constraint and allow x a to vary in the kinematically allowed range x a ,min < x a < 1, 

where 

E p E T e~^ 



2E e E p - E e E T e^ 



(2.40) 



except when we compare to experimental data. There we shall include the correct limits on x a dictated 
by the experimental analysis. From eqs. Q2.37 ) and (2.38), we can express Xb as a function of Ef, T]i, 

and x a : 

x a E e E T e^ 

Xb ~ 2x a E e E p - E p E T e-vi ' ^ l > 
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The inclusive two-jet cross section for ep — ► e + jet x + jet 2 + X is obtained from 

d ^2 d a d =J2 x * F a/e(xa,MZ)x b F b/p (x b ,M{j)^(ab^jets), (2.42) 

^ ' ' a, b 

where 

1 

F a/e (x a ,M a 2 ) = / ^F ah {y a ,M 2 a )F l/e (^) (2.43) 

J Va \y a J 

defines the parton content in the electron, da/dt stands for the differential cross section of the 
partonic subprocess ab — ► jets. The invariants of this process are s = (p a + Pb) 2 -, t = (p a — Pi) 2 , and 
u = (p a — P2) 2 ■ They can be expressed by the final state variables Et, 771, and 772 and the initial state 
momentum fractions x a and x b : 

(2.44) 
(2.45) 
(2.46) 

For the inclusive one-jet cross section, we must integrate over one of the rapidities in ( |2.42| ). We 



s 


QX a X b il/ e rjp, 




t 


= -2x a E e E T e^ = 


-2x b E p E T e- V2 


u 


= -2x a E e E T e^ = 


-2x b E p E T e- m 



integrate over 772 and transform to the variable x a using (|2.37[) . The result is the cross section for 
ep — > e + jet + X, which depends on Et and 77: 



d 2 a 
dETdr] 



J2 f dx a x a F a/e {x a ,M 2 a )x b F b/p {x b , M b 2 ) *f eE l ^-{ab -> jet). (2.47) 

tt J x a ,mm 2x a E e - E T e V dt 



Here, x b is given by (2.41) with r\\ = 77. 



3 Leading Order Cross Sections 

In this section, we consider the leading order contributions to the photoproduction cross section dcr/dt. 
The leading order direct hard scattering process is of 0(aa s ), and the leading order resolved hard 
scattering process is of 0(a 2 ). Since the latter has to be multiplied with the photon structure function 
of 0(a/a s ), both contributions are of the same order 0(aa s ). To this order, it is necessary to calculate 
the phase space of two-particle final states, the tree-level matrix elements for direct photons 76 — > 12, 
and those for resolved photons ab — > 12. The two particles produced in the hard scattering then 
correspond directly to the two jets that can at most be observed. Either particle can be in jet one or 
jet two, so that the cross sections da/dt have to be considered in their given form and with (t <-> u) to 
give the complete cross section da/dt(s,t,u)+da/dt(s,u,t). Furthermore, one has to add symmetry 
factors of 1/n! for n identical particles in the final state. As we will later go from four to d = 4 — 2e 
dimensions to regularize the singularities showing up in next-to-leading order, we already give the 
results in this section in d dimensions. If one is only interested in leading order results, one can always 
set e to zero in this section. 

Section 3.1 contains the phase space for 2 — ► 2 scattering in d dimensions. It is calculated in the 
center-of-mass system of the incoming or equivalently outgoing particles. The matrix elements for 
direct photons are presented in section 3.2. Only one master diagram for the so-called QCD Compton 
scattering process jq — > gq contributes. The gluon initiated process is obtained by crossing. For 
resolved photons, we have four parton-parton master diagrams. The corresponding matrix elements 
are given in section 3.3. In addition, we show a table from which the matrix elements for all other 
processes can be deduced. Section 3.4 contains the direct 77 Born matrix element. 
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3.1 Phase Space for Two-Particle Final States 

We give here a brief sketch of the phase space calculation for the scattering of two initial into two final 



particles in d = 4 — 2e dimensions. We start from the general expression [48] 



d ps<*> = /<■*? n *g!$* L +n -± Pj \ 



(3.1) 



The dijet and single-jet cross sections in eqs. ( |2.42| ) and ( |2.47| ) require partonic cross sections differential 
in the Mandelstam variable t. We therefore insert an additional ^-function with respect to t 



dPS^) 
dt 



/( 2 -) d A ^0r sd U + w " I>) *(* + 2p-pi) 



(3.2) 



before we integrate over the <i-dimensional <5-function leading to 

Next, we choose the center-of-mass system of the incoming partons a, b as shown in figure @, where 




Figure 6: Center-of-mass system for the scattering of two initial into two final partons. 

the outgoing partons are of the same energy and can be described by a single polar angle 6 between 
p\ and p a . We integrate over the azimuthal angle 4> 



dPS< 2 > 
dt 



dExEf-* 'dcos 0(1 - cos 2 9) — 

2 d - 2 7r^r to 



■<J(s - 2>/l£i)<$(t + >/a^i(l - cos6»)) 



(3.4) 



and over the remaining 5- functions. The final result is [149, 50, 51] 



dPS^) 
dt 



1 



/47TSV 



1 



T(l -e) \tu J 8vrs 



(3.5) 
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3.2 Born Matrix Elements for Direct Photons 

For direct photoproduction of two partons, there is only one generic Feynman diagram "fq — ► gq as 
displayed in figure f?|. Unlike deep inelastic scattering, the leading order process already includes the 



Tnnnr 



Figure 7: 2 — ► 2 Feynman diagram for direct photoproduction. 

production of a hard gluon, which balances the transverse momentum of the scattered quark. This 
can either happen in the initial state (left diagram in figure |8|) or in the final state (right diagram in 
figure H) of this so-called "QCD Compton Scattering" . Both diagrams have to be added and squared 



q 7 



innnr 



g q 



Figure 8: Born diagrams for direct photoproduction. 
to give the full matrix element squared. In d dimensions, the result is 

\M\ 2 q ^ gq (s,t,u) = e 2 e 2 q g 2 n i£ T iq ^ gq (s,t,u), 
where the renormalization scale /x keeps the couplings dimensionless and where 



T jq ^ gq (s,t,u) = 8N C C F {1 - e) 



c 1 -*)(----) 

V s uj 



+ 2e 



(3.6) 



(3.7) 



Note that the DIS interference term of the two diagrams is missing as it is proportional to Q 2 ~ 0. 
As discussed in section 2.2, the proton not only consists of quarks but also of gluons giving rise to a 
"Boson Gluon Fusion" process jg — ► qq not shown here. The corresponding matrix element is obtained 
from eq. ( |3.6D by simply crossing the initial quark with the final gluon, or (s <-► £), and multiplying 
by (—1) for crossing a fermion line. Similarly, the contribution for incoming anti-quarks "fq~ — > gq can 
be calculated by crossing (s «-> u). The LO direct matrix elements have been known for quite some 
time m, 52, 53 1 and are summarized in table |2|. 



Process 



Matrix Element \M\ 2 



[\M\ 2 q ^ aq (s,t,uW¥m 



jq^ gq 



¥ 



7<?^ gq 



[\M\ 



ii^m 



[u,t,s)}/[ANc] 



lM^ Q ^ qq ( t, S ,u)]/[8(l-e)N c C F ] 



ia^ qq 



Table 2: Summary of 2 — v 2 squared matrix elements for direct photoproduction. 
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3.3 Born Matrix Elements for Resolved Photons 

For resolved photoproduction of two partons, we have to calculate the four generic parton-parton 
scattering processes qq' — > qq', qq — > qq, qq — > gg, and gg — ► gg. The corresponding Feynman 
diagrams are displayed in figure 0. The first process for the scattering of two non-identical quarks 



a) 



q q 



q q 



-9 9 



b) 



Tnnnr "innny Tnnnr 



c) d) 



-9 



Figure 9: 2 — ► 2 Feynman diagrams for resolved photoproduction. 

q and </ has already been calculated in 1978 by Cutler and Sivers [p4| , shortly prior to the other 
processes p5[ , and was compared to inclusive jet and hadron production at ISR energies. 

Figure p!o|a) shows that for qq' — v (?(/ only the one-gluon exchange between the different quarks in 
the ^-channel contributes. The matrix element squared is given by 



\M\ 2 qq ,^ qq ,{s,t,u) 



with 



fv 4£ T qq i-, qq '(s,t,u) 



s 2 + u 2 



T qq t^ qq ,(s,t,u) = AN C C F I — -2 e 



(3.8) 
(3.9) 



We will split off the coupling constants for the other Born matrix elements \Ai\ 2 as in eq. (| 

Equal quark flavors cannot be distinguished in the final state so that both diagrams in figure |H 
have to be added. In addition to the squares of the individual diagrams already present in the case of 
different quark flavors, the interference term 



T, 



qq^qq\S> *) ^7 



-8C P (l-e) 



ut 



+ e 



(3.10) 



also contributes to the process qq — > qq. 

Final gluons are also undistinguishable and can even couple in a non-abelian way as shown in 
figure |l0|c) . The resulting matrix element for the process qq — ► gg is 



L qq^gg 



(s,t,u) =4Cp(l -e) 



(2N C C F 2AT2 

y ut 



c 



(t 2 + u 



2 es 2 ) 



The diagrams for the process gg — > gg in figure Kd) give 



T gg ^ gg (s, t, u) = S2N£C F {1 - e) 2 3 



ut us 



st\ 



(3.11) 



(3.12) 



All other diagrams can be obtained from the above by simple crossing relations. The complete 



result in d dimensions can be found in [56] and is summarized in table 
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* nnnnr 



g q 



* — (nnnr 



g q 




g g 



g g 




b) 




Figure 10: Born diagrams for resolved photoproduction. 



Process 


Matrix Element |.M| 2 


qq' — ► qq' 


O^C^'M^M 4 ^ 


qq' -► qrg' 


[l^|^^(«,t,5)]/[4iVg] 


<?g -► <?V 


[|M|^_^(M,«)]/[4iV£] 


qq^qq 


[l*C- 


^(5,t,«) + |M|^ 00 /(8,u,t) + |A^|^ M (s,t,«)]/[4iV^/2! 


qq^qq 


[l*C 


_>„,(«, t, a) + |M|^,^ w ,(n, a, t) + |A4|^ M (u, t, s)]/[4iV£] 


qq^ gg 


[|-M|^ M (*,*,«)]/[4iVS]/2! 


qg -► w 


-0^|^ w («,a,«)]/[8(l-e)JV^C F ] 


95 "^95 


-[|^|^ M (*,«,s)]/[8(l-e)^C F ] 


55^ 99 


[|^|^ aa (,,t,n)]/[16(l- £ )^^] 


55 ^55 


[|^|^ Q ( S ,t,u)]/[16(l- £ )^^]/2! 



Table 3: Summary of 2 — > 2 squared matrix elements for resolved photoproduction. 
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3.4 Born Matrix Element for Direct 77 Scattering 

Since the gluon has no electromagnetic charge, direct 77 scattering can only proceed through the 



process 77 — ► qq as shown in figure 11. The corresponding matrix element is given by 



Figure 11: Born diagram for direct 77 scattering. 



\M\ 2 ^ qq (s,t,u) = e 4 eg/i 4e T 77 ^gg(s,t,u), 



where 



<■ 77-xjg 



(s,t,u) = 8N c (l-e) 



(1 



u t 

t u 



2e 



(3.13) 



(3.14) 



It can be obtained from the direct photoproduction matrix element for ^q — ► gq in |3.6| through crossing 
of s — > t, multiplication with (-1) for the crossing of a fermion line, and replacing the strong coupling 
constant g 2 by e 2 e 2 and the color factor Cp by 1. 



4 Next- To-Leading Order Cross Sections 

As in any calculation of next-to-leading order cross sections, we have to calculate two types of cor- 
rections in photoproduction: virtual and real corrections. For the direct case, these corrections are of 
0(aa 2 ), and for the resolved case, they are of 0(o? s ). Throughout this section, direct and resolved 
photon results will be presented separately, but in a completely parallel way. Whereas the direct 
contributions have already been published |37], g|, the resolved contributions are presented here for 
the first time. 

The virtual corrections consist of the interference terms of one-loop graphs with the Born graphs 
calculated in section 3 and will be given in section 4.1. We can still use the same phase space for two 
final partons as in the last section. The direct matrix elements are contained in section 4.1.1, those for 
resolved photons in section 4.1.2. Again, we display only the relevant master diagrams, from which 
all subprocesses can be deduced through crossing. Section 4.1.3 contains the virtual corrections for 
direct 77 scattering. 

The real corrections are derived from the integration of diagrams with a third parton in the final 
state over regions of phase space, where this third parton causes singularities in the matrix elements. 
Different methods can be employed here. We choose the phase space slicing method with an invariant 



mass cut [44] in section 4.2. Alternatively, the subtraction method could be used |6^, |58|, y]. We 
calculate the three particle phase space in two different versions for final state singularities (section 
4.2.1) and initial state singularities (section 4.2.4). These sections are followed by the calculation of 
the matrix elements for final state (sections 4.2.2 through 4.2.3) and initial state corrections (sections 
4.2.5 through 4.2.8). The real corrections for direct 77 scattering are presented in section 4.2.9. 

Finally, we demonstrate the cancellation of the infrared and collinear singularities, that show up 
in the virtual and real corrections separately, in section 4.3. The ultraviolet singularities in the virtual 
corrections are removed by counter terms in the renormalization procedure. Remaining collinear 
singularities in the initial states can be absorbed into the photon and proton structure functions. 
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4.1 Virtual Corrections 

The central blobs in figures p] and U can also contain loop corrections to the leading order matrix 
elements in figures || and |l0[ Up to 0{aa 2 ), the interference terms between the direct Born diagrams 
T and the direct one-loop diagrams V have to be taken into account according to 



\M\ 2 (0{aa 2 s )) = \^fao^T + ^ aalV\ 2 

= aa s T + 2aa 2 s W*T + 0{aa 3 s ). (4.1) 

Similarly, the interference terms between the resolved Born diagrams T and the resolved one-loop 
diagrams V have to be taken into account up to 0{p? s ). 

The inner loop momenta are unconstrained from the outer 2 — ► 2 scattering process. We can 
therefore still make use of the phase space calculated in section 3.1, but have to integrate over the 
inner degrees of freedom. At the lower and upper integration bounds, we encounter infrared (IR) and 
ultraviolet (UV) divergencies. Using the dimensional regularization scheme of t'Hooft and Veltman 
[p9J ], we integrate in d = 4 — 2e dimensions, thus rendering the integrals finite and keeping the 
theory Lorentz and gauge invariant. The IR and UV poles show up as terms oc 1/e 2 and 1/e. As 
a consequence, the complete higher order cross sections have to be calculated in d dimensions up to 
0(e 2 ) to ensure that no finite terms are missing. 

The ultraviolet divergencies coming from infinite loop momenta can be removed by renormalizing 
the fields, couplings, gauge parameters, and masses in the Lagrangian, since QCD is a renormalizable 
field theory. This is done by multiplying the divergent parameters in the Lagrangian with renormal- 
ization constants Z% and expanding up to the required order in the coupling constant g. The resulting 
counter terms then render the physical Green's functions finite. In addition to the 1/e poles, universal 
finite contributions 

- - 7B + ln(4vr) (4.2) 

are included into the counter terms according to the modified minimal subtraction (MS) scheme IJBCR , 
which we employ here. The coupling constant g is kept dimensionless by multiplying it with the 
renormalization scale \x 

9 - gif- (4.3) 

If eq. ( |4.3[ ) is expanded in powers of e and combined with single poles of the type in eq. ( |4.2[ ), the 
renormalization procedure will produce an explicit logarithmic dependence on the renormalization 
scale fj, 

l(*£)\l +hL *£, (4.4) 

e y s J e s 

that cancels the first term of the expansion of the running coupling constant 
. 2 , 12tt / 6(153 -19i\Mln(ln£)\ 

"-^- (s.-awg i 1 - (33 W ln g J*"" < 4 ' 5 > 

a s is given here in two-loop approximation as appropriate in NLO QCD, and Nf is the number of 
flavors in the qq loops. We choose ^ = 0{Et) as for the factorization scales. 

Infrared divergencies arise from small loop momenta in the virtual corrections. They have to cancel 
eventually against the divergencies from the emission of real soft and collinear partons according to 
the Kinoshita-Lee-Nauenberg theorem |pl| . Finally, one obtains finite physical results in the limit 
d^A. 
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Figure 12: Virtual diagrams for direct photoproduction. 
gluon, and ghost loop. 



The circle in diagram a) denotes a quark, 



4.1.1 Virtual Corrections for Direct Photons 

The one-loop corrections to the left diagram of figure || are shown in figure [jl| and can be classified into 
a) self-energy diagrams, b) propagator corrections and box diagrams, and c) vertex corrections. They 
contain an additional virtual gluon, which leads to an extra factor a s . Similar diagrams are obtained 
for the right diagram of figure [8|. As in leading order, the diagrams for the process 7<? — ► qq can be 
obtained from figure [TJ by crossing the initial quark and the final gluon or equivalently (s «-► t) and 
multiplying by (—1) for the crossing of a fermion line. However, the virtual corrections can in general 
contain logarithms ln(x/s — irf) and dilogarithms Li2(x/s), where x denotes different Mandelstam 
variables before and after crossing. The ^-prescription for Feynman propagators then takes care of 
additional terms of tt 2 , that arise in the quadratic logarithms with negative argument. 

The virtual contributions have been well known for many years from e + e~ — > qqg higher order QCD 
calculations [62, 63]. For the corresponding photoproduction cross section, one substitutes Q 2 = 



and performs the necessary crossings. The result can be found in [49, |33|, 57 1. We have also compared 
with the results in [J5l| for deep inelastic scattering eq — ► e'gq and eg — ► e'qq, which can be expressed 
by the invariants s, t, and u after setting Q 2 = 0. For completeness and for later use, we write the 
final result in the form 
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The expression V iq ^ gq (s,t,u) is given by 
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The contributions for incoming anti-quarks and gluons ~fq — ► gq and 75 — ► gg could be obtained 
according to table 0, if the imaginary parts were included above. The result for anti-quarks turns out 



to be identical to eq. (4.7), and for the gluon initiated process one obtains 
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All UV poles have been canceled by counter terms through the renormalization procedure, and the 
remaining IR poles are contained in the first three lines of eqs. ( f4.7[ ) and ( |4.8[ ). Contrary to some of the 
finite terms, they are proportional to the LO Born matrix elements T^ q ^ gq and T ig ^ q q. The explicit 
dependence of the virtual corrections on the renormalization scale fi is contained in the logarithms 
proportional to the Born matrix elements and Nq and Nf, respectively. All terms of 0{e) and higher 
have been omitted since they do not contribute in the physical limit d — ► 4. 



4.1.2 Virtual Corrections for Resolved Photons 

For the resolved virtual corrections, we show in figure ^ only the one-loop diagrams for qq' — > qq' and 
classify them again into a) self-energy diagrams, b) propagator corrections and box diagrams, and c) 
vertex corrections. Due to the additional gluon, these contributions are again one order higher in a s 



than the Born terms as in the direct case. The complete set of diagrams can be found in [56] as well 
as the results 
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where Q 2 denotes now an arbitrary scale. For the diagrams in figure ^, we obtain 



V qq >^ qq '(s,t,u) 
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Figure 13: Virtual diagrams for resolved photoproduction. Only the one-loop corrections to qq' — > qq' 
are shown. The circle in diagram b) denotes a quark, gluon, and ghost loop. 
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The contributions for incoming quarks and anti-quarks of different flavor can be obtained from table 
H as well as the contribution for identical quark flavors. There one also has to include the interference 
term 
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For incoming antiquarks, the Mandelstam variables (a *-* u) have to be crossed. The first lines 
proportional to the color factors Cf and Nc contain the IR singular terms in the two processes 
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discussed above. Furthermore, the complete virtual corrections are proportional to the Born matrix 
elements. This is not true for the case qq — ► gg, where the virtual corrections are given by 
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Here, only parts of the IR singular terms are proportional to the Born matrix element T c 
finite contributions have been put into the function 
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Finally, the process 5(7 — ► 55 gives 
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where, once more, not all IR singularities are proportional to the complete Born matrix element and 
the finite terms are given by the function 



f d (s,t,u) = N c 



2{t 2 + u 2 )\ , 2 . , fis(t 3 + u 3 ) \ ,..,, . 



+ 



' Atu Ut 2 + u 2 
3^ ~~ 



tu 



t 2 V? 

ft 2 u r 



(4.15) 



l(s) 



1-7T' 



+ - 



iVf 



t 2 +M 2 



£u 



'l0t 2 + -u 2 16 tu 
3 tu 3s 2 

Z(i)Z(u) + 2 - vr 2 



2 W( s ) 



S 2 + tU 2 



t'U 



/ 2 (.) 



All UV singularities have already been absorbed into the Lagrangian, and terms of 0(e) and higher 
have been neglected. 

As mentioned previously, the crossing of Mandelstam variables according to table |3| may change 
the signs of the arguments in the logarithms. This is accounted for in the definition of 
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where x may be any of the variables s, t, or u and is normalized to Q 2 = max(s,t,u). Due to the 
^-prescription in the propagators, terms bilinear in the logarithms may give an extra n 2 : 
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4.1.3 Virtual Corrections for Direct 77 Scattering 

The virtual diagrams contributing to direct 77 scattering are shown in figure |lj. As has already been 
seen in the Born case, direct 77 — > qq scattering is intimately related to jg —* qq scattering. Again, 
we have to replace the strong coupling constant by its electromagnetic counterpart resulting in 



The expression V yi ^ g q(s, t, u) is given by 
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Note that due to the abelian structure of QED, the color class Nc, which is present in eq. (4.8) and 
arises there from the triple-gluon vertex, does not appear here. 
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Figure 14: Virtual diagrams for direct 77 scattering. 



4.2 Real Corrections 

For the calculation of the hard scattering cross section in next-to-leading order, we must include all 
diagrams with an additional parton in the final state. The four-vectors of these subprocesses will be 
labeled by p a Pb —> P1P2P3 , where p a is the momentum of the incoming photon or parton in the photon 
and ph is the momentum of the incoming parton in the proton. The invariants will be denoted by 
s ij = (Pi +Pj) 2 an d the previously defined Mandelstam variables s, t, and u. For massless partons, the 
2 — ► 3 contributions can contain singularities at sy = 0. They can be extracted with the dimensional 
regularization method and canceled against those associated with the one-loop contributions. 

In order to achieve this, we go through the following steps. First, we calculate the phase space 
for 2^3 scattering in d dimensions and factorize it into a regular part corresponding to 2 — ► 2 
scattering and a singular part corresponding to the unresolved two-parton subsystem. Next, the 
matrix elements for the 2^3 subprocesses are calculated in d dimensions as well. They are squared 
and averaged/summed over initial/final state spins and colors and can be used for incoming quarks, 
antiquarks, and gluons with the help of crossing. 

One can distinguish three classes of singularities in photoproduction (Q 2 ~ 0) depending on which 
propagators in the squared matrix elements vanish. Examples for these three classes are shown in 
figure [l5|. The X marks the propagator leading to the divergence. In the first graph, it is the invariant 
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Figure 15: Three-body diagrams with a) final state, b) photon initial state and c) proton initial state 
singularities. 

S23 given by momenta of the final state, which causes the singularity. Therefore, this divergence will 
be called final state singularity. The second graph becomes singular for —t a \ = s a \ = (p a + Pi) 2 = 0, 
when the photon and the final quark momentum are parallel. This is the class of photon initial state 
singularities. In the third graph, the singularity occurs at — t^ = Sb3 = (pb +P3) 2 , where pb is the 
initial parton momentum in the proton. This stands for proton initial state singularity. Since resolved 
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photons behave like hadrons, they produce similar initial state divergencies. The first class is familiar 



from calculations for jet production in e + e annihilation [63], the third class from jet production in 
deep inelastic ep scattering (Q 2 ^ 0) |J5l| . The second class occurs only for direct photoproduction 



When squaring the sum of all relevant 2^3 matrix elements, we encounter terms, where more than 
one of the invariants become singular, e.g. when one of the gluon momenta p 3 — > so that S23 = and 
—tb3 = Sfc3 = 0. These infrared singularities are disentangled by a partial fractioning decomposition, 
so that every term has only one vanishing denominator. This also allows the separation of the different 



classes of singularities in figure 15, 



It turns out that the results for direct photoproduction are always proportional to the LO cross 
sections involved in the hard scatter 

F,I,J lb ^ 123 = KT lb ^ 12 , (4.21) 

where F, I, and J denote final state, photon initial state and proton initial state contributions. For 
resolved photoproduction, this is no more true in general but only for the quark-quark scattering 
subprocesses with a less complicated color structure than those involving real gluons. For gluonic 
processes, only parts of the leading order cross sections can be factorized 

F,I,J ab _> 123 = J2K^T^ 12 , (4.22) 

i 

where 

T^i2 = J2 T !Li2 (4-23) 

i 

is again the full Born matrix element. 

As the last step, the decomposed matrix elements have to be integrated up to Sy < ys. y char- 
acterizes the region, where the two partons i and j cannot be resolved. Then, the singular kernels K 
and K W produce terms oc 1/e 2 and 1/e, which will cancel against those in the virtual diagrams or be 
absorbed into structure functions, and finite corrections proportional to In y, Iny, and y°. For small 
values of y, terms of 0(y) can be neglected. In the following, we shall give the results for the different 
classes of singularities separately. 

4.2.1 Phase Space for Three-Particle Final States 

For the real corrections, we have to consider all subprocesses, which have an additional third parton 
in the final state attached to the 2 — ► 2 scattering process. We therefore calculate the phase space 
for 2 — ► 3 scattering in this section and choose as the coordinate system the center-of-mass system 



of partons 1 and 3 as shown in figure 16. The angle between p a and p\ is called 9, the azimuthal 
angle between the planes defined by p a and p\ and p a and p 2 is called 0, and the angle x* between 
p a and p 2 is defined in the overall center-of-mass system of the incoming partons a and b denoted 
by an asterisk (*). As we have to accommodate a third final particle in the scattering process, the 
Mandelstam variables in this section 

(4.24) 

2 PlP3 , (4.25) 

u - (p a -P2)' z (4.26) 

differ slightly from the previously used ones, but still satisfy the relation s + t + u = for massless 
partons. In the limit of soft (p 3 = 0) or collinear {p 3 || pi) particle emission, they can however be 
written in a form similar to 2 — ► 2 scattering as t = (p a —p-\) 2 and u = (p a — p 2 ) 2 - Here Pi = Pi + P3 
represents the four-momentum of the recombined jet, and p 2 is the four-momentum of the second jet. 
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Figure 16: Center- of-mass system of partons 1 and 3 defining the kinematic variables for the scattering 
of two initial into three final partons. x* * s defined in the overall center- of-mass system of partons a 
and b. 



We start again from the general expression 48 
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and describe the final state singularity with the new variable 
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We can then insert two additional (5-functions with respect to t and z' 
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before we integrate over the 5(pf) and the space-like components of the d-dimensional 5-function to 
eliminate p 3 . The resulting expression is 
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We now decompose pi into its energy E\ and angular components 6 and in the center-of-mass system 
of partons 1 and 3, and p2 into its components E%, x*, and §\ in the overall center-of-mass system 
with the result 
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Integrating over the remaining 5- functions and the trivial azimuthal angle $\ up to 2ir, we arrive at 

(4.32) 
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Finally, we can factorize this three particle phase space into 

dPS (3) = dPS (2) dPS (r) , 
where 
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is the phase space for the two observed jets with momenta pi and p2 already calculated in section 3.1 
and 

' : dfx F (4.38) 



dPS^ 



4^\ £ r(l-e) 



s J T(l - 2e) 16vr 2 1 - 2e 
is the phase space of the unresolved subsystem of partons 1 and 3. The integration measure is 
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The full range of integration in dPS^ is given by z' £ [0, — t/s], b € [0, 1], and 4> £ [0, 7r]. The singular 
region is defined by the requirement that partons p\ and P3 are recombined, i.e. p^ = or j?3 parallel 
to pi, so that Si3 = (pi +P3) 2 = 0. We integrate over this region up to S13 < ys, which restricts the 
range of z' to < z' < min{— t/s, y} = i/f- 

4.2.2 Final State Corrections for Direct Photons 

The real corrections to the QCD Compton Scattering process of figure [7| arise from two different 
mechanisms. Either an additional gluon can be radiated from the quark or the gluon leading to the 



left diagram in figure 17, or a quark- antiquark pair is emitted by a gluon as shown in the right diagram 






9 7 



q q 



q,q' 



Figure 17: 2^3 Feynman diagrams for direct photoproduction. 
of figure [l?]. Both cases lead to an extra factor of a s , when the matrix elements are squared, so that 
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Process 



Matrix Element \M\ : 



JM^ 



ii -> qgg 



•jm_ 



[8,t,u)]/[4Nc]/2\ 



1Q -> 955' 



[l-MI 



■yq-^qgg 



(s,t,u)}/[AN c }/2\ 



\\M^ qqq -(s,t,u)]/[W c ] 



7g -> ggg 



7ff -" 599 



;M£ 



( S ,t,«)]/[8(l-e)^ c C F ] 



Table 4: Summary of 2 — ► 3 squared matrix elements for direct photoproduction. 



all real corrections discussed in the following are of O(ao%). The outgoing quark-antiquark pair can 
be of the same or different flavor than the incoming quark. The diagrams for incoming anti-quarks or 



gluons can be obtained from figure 17 by crossing a final quark or gluon line with the incoming quark 
line. The corresponding matrix elements can then be obtained from table ||. 

All possible topologies and orders of outgoing particles have to be considered for the full matrix 
elements of the processes, although they are not shown here explicitly. The complete result in d 
dimensions can be found in [Ej| in a very compact notation not suitable to isolate the singularities. 



We therefore re-calculate the matrix elements with the help of REDUCE [64], check their sums with 
the results in Eg], and keep only those that have singularities in a final state invariant sy = (pi +pj) 2 . 
This can either be a soft or collinear gluon leading to a quadratic pole (left diagram of figure |17|) or 
a collinear quark-antiquark pair leading only to a single pole (right diagram of figure 17). With the 
help of the singular invariant 

z> = ™ , (4.40) 

PaPb 

we can approximate the nine other invariants describing the 2 — ► 3 scattering process and express 
them through the 2 — ► 2 Mandelstam variables s, t, and u and the variable b = 1/2(1 — cos 9) as 
defined in section 4.2.1: 

(4.41) 

(4.42) 

(4.43) 

(4.44) 

(4.45) 

(4.46) 

(4.47) 

(4.48) 

(4.49) 

(4.50) 





S 


PaPb = 


2 




s ~ t i. 


PaPl = 


- 2~ b 




s —u 


PaP2 = 


27 




8-t, 


PaPZ = 


- IT' 1 - 6 ' 




s —u. 


PbPl = 


= b 

2 s 




s-t 


PbP2 = 


27 




s ~ u ,-. 


PbPZ = 


- 2-'- 1 -"- 




s. 


P\P2 = 


= -b 
2 




s i 


PlP3 = 


= -z 
2 


P2P3 = 


= Ul-b) 



For the subprocess jq — ► qgg, where a gluon becomes soft or collinear to the other outgoing gluon or 
the quark, the approximated matrix element takes the form 



\M\ 



2 
11^199 



(s,t,u) 



H sz' 



4C F (1 -&)(l-e)- 2 + 



-2t/s 



z' -t/s(l-b) 



(4.51) 
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-AN C 



-t/s 



-u/s 



+ 2-6 + 6' 



, z' - t/s(l -6) z' - u/s(l -6) z' + (1 - b) 

with an identical result for incoming anti-quarks. The four-fermion diagram only produces a divergence 
from the collinear quark- antiquark pair, which can have Nj flavors 

\M\*^ n9 (8,t,u) = e 2 e yn 6£ ^N f [l - 26(1 - 6)(1 + e)]T jq ^ gq (s,t,u), (4.52) 

and the gluon-initiated process 75 — ► gqq can have a gluon that becomes soft or collinear to the quark 
or anti-quark 



\M\ 2 g ^ gqq (s,t,u) 



e 2 ey^-, 

H sz 



-2N C 



4C F ((l-b)(l-e)-2 
-t/s 



2 



z> + (1 - 6) 

—u/s 



z' + (l-b) z'-t/s(l-b) z'-u/s(l-b) 



T- 



-yg^qq 



(4.53) 



(s,t,u). 



Finally, we integrate the above matrix elements over the singular phase space dPS^ given in 
section 4.2.1. The necessary integrals are given in appendix A. The integration over in dPS*- 7 *' is 
trivial, as the matrix elements do not depend on this variable. This is true for all direct and resolved 
photoproduction processes and for all final and initial state contributions, contrary to the case of 
deep-inelastic scattering. The result is 



dPSM|M|^ 123 (M,u) = e^sV'g 



2 2 2 4e a s f^^ 2 \ r (!-e) 



r(i - 2s) 



F 7 b-*l23(s,t,ll). 



(4.54) 



The functions .F 7 b-»i23 are 


given by 




^jq— *qgg[ s i *> u ) = 


~ 2CF {h- 


2e \ s 




vr 2 7 
3 + 2 


1 9 -t 

- - In 2 h 2 In y 

2 s y 



2Li 2 f — — j - In 2 y F - ~ In y F 



3 

2 



„ / 2 1 (11 . -t . -u 

1, ,-t , , -t T /?/fs\ 67 11, , o 

--In 2 — +2 In y F In 2Li 2 ^- - — + — lny F + In 2 y F 

2S S !<;<>•> 



V * / 



9 





+ ln 2 ^-Iln 2 ^ + ^+2Li 2 f^V 
— u 2 s 3 \ u / 1 


r yq^qqq\S,t,U) — J\ f 


r 1,1, 5" 

1 — In «f 

3e 3 y 9 


l'yq-^gq\ s i ^1 u )j 


ig^gqq\ s i *j w / == 


C B 


/ 2 3 2vr 2 , 
\e z e 3 


7 -2ln 2 y F -3lny F \ 



-iq^gq 



(s,t,u), 



(4.55) 
(4.56) 



N c (1 



2 Vfs 



2 2/fs 



- In V In + ln z ^— + In" 4 — 2 ln z y F - - ln z — 

2 \e \ s s J —t —u 2 s 



Z±4.mAJvz£\+*iAJm 



.Il n 2_ + 2 Li 2 (^l)+2Li 2 
2 s V t J \ u 



yg^qq \ s > ^1 u ) 



(4.57) 



Terms of 0{e) and of 0{y F ) have been neglected. The Spence-functions or dilogarithms Li 2 resulting 
from the calculation are written down explicitly even though they are also of 0(y F ). The complete 
approximated 2 — > 3 matrix elements are proportional to the corresponding Born matrix elements 
T 7 (,^i 2 . This is also true for the exact matrix elements in four dimensions as calculated by Berends 
et al. [|65|1, but not for the exact d-dimensional matrix elements 
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4.2.3 Final State Corrections for Resolved Photons 

For the four different generic Born diagrams of resolved photoproduction in figure 0, the real corrections 
all arise from an additional gluon in the final state as shown in figure 111. They are of C(a^). All 




a) 



-9 9 



b) 



Tnnnr "innny Tnnnr 



c) 



Tnnnr 



9 
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d) 



-9 



9 
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Figure 18: 2 — ► 3 Feynman diagrams for resolved photoproduction. 

other diagrams can be obtained from figure ITS by crossing one or two outgoing gluon lines into the 
initial state, which leads to the gluon initiated processes with quarks in the final state, or by crossing 
an incoming and and outgoing quark line, which leads to the processes with incoming anti-quarks. 
Process c) is symmetric under the interchange of the three gluons and under the interchange of the 
quark and anti-quark. Process d) is completely symmetric under the interchange of any of the five 
gluons. The complete list of matrix elements is written down in table || 



Process 


Matrix Element \M\ 2 


qq' -> qq'g 


[\M\ 2 q ^ qqlQ (s,t,u)]/[AN 2 ] 


qq' -> qq'g 


[\M\ 2 q ^ qqlq {u,t,s)]/[AN^\ 


qq -> q'q'g 


[\M\^ qq ,{t,s,u) + \M\ 2 q ^ qqq (s,t,u)}/[W£,} 


qg -> qq'q' 


-[\M\^ qqq (t,s,u)]/[8(l-s)N^C F ] 


qq -> qqg 


l\M\' 2 qql _> qq , q (s,t,u) + \M\* q r_» qqlq (s,u,t) + \M\ 2 q ^ qqq (s,t,u)]/[4N 2 }/2\ 


QQ -► QQ9 


[\M\ 2 qq ,^ qqlq (u,t,s) + \M\ 2 q ,^ qqlq (u,s,t) + \M\ 2 qq ^ qgq (u,t,s) + \M\ 2 qq ^ qqq (s,t,u)]/[4N 2 ] 


Q9 ~> QQQ 


-MM\ 2 qq qq {t, s, u)]/[8(l - e)N 2 C F ]/2\ 


QQ ~> 999 


[\M\ 2 qq ^ qqq (s,t,u)]/[AN 2 ]/3\ 


Q9 -> Q99 


[-\M\ 2 q ^ qqq (t, S , U )/3 + \M\ 2 q ^ qqq (s,t,u)]/[8(l-e)N 2 C F ]/2\ 


Q9 -> Q99 


l-\M\ 2 qq -_ qqq (t, u, s)/3 + \M\ 2 qq ^ qqq (u, t, S )]/[8(1 - e)iV£0]/2! 


99 -> QQ9 


[-NI^ O00 (i,^«) + l^|^ 0O - (5,t,«)]/[16(l- £ )^C|] 


99 -► 999 


[|^|^ M ( a ,t,«)]/[16(l- e )^crj,]/3! 



Table 5: Summary of 2 — ► 3 squared matrix elements for resolved photoproduction. 
In four dimensions, these matrix elements were first calculated by Gottschalk and Sivers [66] and 



by Kunszt and Pietarinen [67]. The result in d dimensions was given in a compact form by Ellis and 
Sexton [56]. We calculate the matrix elements again with the help of REDUCE pi], check with the 
result in [56], and keep only terms singular in the final state invariant z' . The result can be expressed 
through this variable z', b, and the Mandelstam variables s, t, and u. For quark-quark scattering of 
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different flavors in diagram [lqa) , we can have a soft or collinear gluon in the final state 

-2 —t/s —u/s 



\M\ 



qq' ^qq' g\ S i *) U ) 



sz 



(2C F - N c ) 



+ 



+ 



+C F (l-6)(l-e)-2 + - 



2' + (1-6) z'-t/s(l-b) z'-u/s(l-b) 
-2u/s 



z' - u/s(l - 6) 
For equal flavors, there is additionally the interference contribution 



'-qq'^qq 



\M\ 



2 

qq^qqg 



(s,t,u) 



ffV 6£ — 

sz 



(2C F - N C ) 



+ 



-t/s 



(s,t,u). 



+ 



(4.58) 



-u/s 



+C F [(l-b)(l-e)-2 + - 



z' + (1 - b) z' - t/s(l -b) z> - u/s(l - b) 
2 



z' + (l- b) 



- qq-^qq 



s,t,u) 



(4.59) 



as the final quark lines in diagram p"g|p) are undistinguishable and can be interchanged. The crossed 
diagram of quark-antiquark annihilation has an additional singularity, if the final quark-antiquark pair 
with Nj possible flavors becomes collinear 



\M\ 



1 



qq^qqg 



s,t,u) = gy b£ —N f [1 - 26(1 - b)(l + e)]T qq ^ gg (s,t,u). 



SZ' 



(4.60) 



The complete matrix elements given above are proportional to the corresponding Born matrix elements 
T a b^i2- Quarks and anti-quarks can, of course, also annihilate into three final gluons as shown in figure 
|lqc), where any pair of two gluons can produce a soft or collinear singularity: 



\M\ 2 qq ^ ggg (s,t,u) 



sV 6e 



12 



sz' 
+4N C C F 
I —u/s 



-N c (b 2 -b + 2)T qq _ gg (s,t,u 



(4.61) 



+ 



z' - u/s{\ - 6) 
1 



1-e 



' t 2 + u 2 -es r 



tu 



+ 



z' + (l- 6) 
-t/s 



(1-e) 



tu 

t 2 +U 2 - £S r 



2N C C F - 2N C ^ - N, 



tu 



N, 



z< - t/s(l - 6) 



(1-e) 



' t 2 + U 2 -£S 2 \ ( 

Tu 1' 



c ^ 



2 t 2 + u 2 \ 
WcC F -2N c ^-N c -t 



tu 



t 1 



Here, only parts of the Born matrix element T qq ^ gg can be factorized. The process qg — ► qgg can be 
obtained by crossing a final gluon with the incoming anti-quark or (s *-* t). In this case, either of the 
outgoing gluons can also be radiated from the outgoing quark leading to 



\M\ 2 qg _ qgg (s,t,u) 



ffV 6£ — 
sz 

-AN C C F 



C F ((l-b)(l-e)-2)T„ g (s,t,u) 



-t/s 



z' -t/s{l-b) 



-N 2S +U 



(1-e) 



V C*" 



+ 



t 2 
-u/s 



z' - u/s(l - 6) 
1 



(1-e) 



s 2 + u 2 - et 2 



Sll 



s 2 + u 2 - et 2 



Sll 



z' + {l- 6) 



(1-e) 



' s 2 + u 2 -et r 



Sll 



Sll 



(2N C C F -2N 2 ^ 



su 



(4.62) 



- 2+2 t 



u 



2N C C F -2N C - T -N C - T 






f 2 



t 2 



(2N C C F -2N 2 --N 2 - 
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where again only parts of the Born matrix element T qg ^ qg can be factorized. If also the incoming 
quark is crossed into the final state, the resulting quark-antiquark pair in the final state with Nf 
flavors can produce the collinear singularity 

1 



\M\ 



2 - ( 
gg-*qqg y 



s, t, u) = gy b£ —N f [1 - 26(1 - 6)(1 + e)] T gg _ gg (s, t, 



u 



SZ' 



(4.63) 



proportional to the full Born matrix element. Finally, the five-gluon process in figure [T§|d) can have 
three different pairs of soft or collinear gluons in the final state with the result 



\M\ 2 gg ^ ggg (s,t,u) 



SZ 

+16N^C F 
( -t/ s 



-N c {b 2 -b + 2)T gg ^ gg {s,t,u) 



+ 



z> - t/s(l - b) 

—u/s 
z' - u/s(l - b) 



(l-ef 



2su 8 + u 

+ 



(i- £ yiz- — + 



t 2 s 2 u 2 

2st s 4 + t 4 



?t 2 



+ 



z' + (1 - b) 



{l-ef 



2tu r + u 4 



t 2 u 2 



(4.64) 



Only parts of the leading-order matrix element T gg ^ gg can be factorized. 

Next, the approximated matrix elements have to be integrated over the final state singularity z' in 
phase space up to the invariant mass cut yp. Terms of order 0(yp) can be neglected if yp is chosen 
sufficiently small. We use the phase space for the unobserved subsystem dPS' r ) from section 4.2.1 and 
find 



dPS^\M\ 2 b _ tl23 (s,t,u) = g 4 t i 



4 ie a s ( Ann 2 



r(i-e) 

2tt \ Q 2 J T(l - 2e) 
similar to the direct case. The functions -F a &— >123 are given by 



aft^l23 



(s,t,u) 



(4.65) 



qq'—>qq'g\S> '■j ^) 



■2i2(,-:\r)i-^(-2l(s)+l(f) + l(u))-^ 



-2l z ( — 

VF 



+ 1 ' 2 ( 7t) + l ' (—)) + \{~ 2l \s) + '*(*) + ^)) + 2Li 2 



-Li 2 



VfJ Wf 

VfQ 2 



VfQ 2 



t 



Li 2 



VfQ 2 



+2( V ( i + 1(3 - 2l{u)) + \l 2 {u) + l~l 2 {— 
2e 2 2 \y F 



vr 2 

2Li 2 

3 



2 lnVF 



VfQ 2 



■-qq'-*qq 



(s,t,u), 



(4.66) 



F, 



qq^qqg 



[s,t,u) 



2(2(V-.V ( ,,(-l(-2/( S )+Z(t) + /(«)) -i 



-2/ 2 



IJF 



t 



II 



\ 



1 



-I-/ 2 \— f )+1 2 [— f )) + \{-U 2 {s) + l 2 (t) + l 2 (u)) + 2Li 2 



VfQ 2 



-Li 2 



VfQ 2 



Li 2 



VfQ 2 



1 J 1-2^ , 7 _/2 



+2C F ( 7 2+-(3-2/( S )) + -Z 2 ( S ) + 



s 



In y F 



TV 



2Li 2 



j/fQ 2 



r, 



qq-^qq 



(s,t,u), 



(4.67) 
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rqq^ggg\S-,t,U) 
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11, 5 
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3e 3 9 



. I 2 11 11 , 67 2vr 
3^c ( — + — - y InyF + y 



e 2 ' 3e 



QQ—>99\ S, > i U n 

2- 



1 qq-^gg 
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+6N C C F 
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(4.68) 



+ (~Ks) + 21(a) - 2/ 2 (— ) + / 2 (.s) - ILio 



VfQ 2 



2\2 



iV, 



C" 



(t 2 + u 2 ) 
uts 2 



+ {-h{t) + 2l{t)-2l 2 (—) +/-(/)- ILi 2 



»*7-£ 



„2\ 



+ ( 2i w(<-£) +a( ^(M 



9f0 2 
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^iK-s 
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(4.69) 



+ (-^(s) + 2 ^( s ) " 2/2 (— ) + 1 2 (») - >LJ2 
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\ \ u su \ s su 
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99^999 



[s,t,u) 



\r ( 2 11 11, 67 2tt 2> 
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[s,t,u) 



+48N£C F 
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2/fQ 2 
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3- -r^ + 
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w 2 s 2 




2/fQ 2 
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u 2 t 2 s 2 

1 
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3 ^ + 
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2/fQ 2 
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(4.72) 



All terms of 0(e) have been omitted since they do not contribute in the physical limit d — ► 4. 

As in the virtual corrections, we account for sign-changing arguments in the logarithms before and 
after crossing with the definition of 



l(x) = In 



Q 2 



(4.73) 



where x is any of the Mandelstam variables s, t, and u and Q 2 is an arbitrary scale chosen to be 
Q 2 = max(s,i, u). However, there are no additional terms of tt 2 in the real corrections so that 



Z 2 (x)=ln 2 



Q 2 



(4.74) 



for x > as well as for x < 0. 



4.2.4 Phase Space for Three-Particle Final States Revisited 

Having calculated all real corrections coming from final state singularities, where the singular variable 
was S13 = (pi +pz) 2 , we now turn to the singularities that arise in the photon and proton initial states 
in the variables — t a ^ = s a s = (p a + p^) 2 and —t b 3 = s b 3 = {pb + Pz) 2 - The singular variables are then 
defined as 



PaP3 



PbP3 



and z >» = ™± (4.75) 

PaPb PaPb 

respectively It is convenient to calculate the three-body phase space in the same center-of-mass 
system of the two final state particles now called 1 and 2 as in the case of final state corrections (see 



figure 16), where now p% is defined in the overall center-of-mass system of partons a and b. However, 
we will not choose the angle 9 between the hard jet and the incoming photon or the related variable 
b = 1/2(1 — cos 6) as the second independent variable. Instead it is more convenient to choose the 
fraction of the center-of-mass energy going into the hard subprocess ab — > 12 



P1P2 ^ rv rl , P1P2 _ [Y nl 
G [X a , 1\ and z b = <E [X b , 1J, 



(4.76) 



PaPb PaPb 

respectively In this way, the variables z a and z b describe the momentum of the third unobserved 
particle p$ = (1 — z a )p a and p% = (1 — z )pb- They are bounded from below through the fractions of 
the initial electron and proton energies transferred to the partons with momenta z a p a and z b pb in the 
hard scattering 



X a 



PlP2 

kpb 



and Xu 



P1P2 
PaP ' 



(4.77) 
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As the three-particle phase space can be calculated for photon and proton initial state singularities 
in exactly the same manner, we only consider the case of the photon initial state in the following. The 
Mandelstam variables differ slightly from those used in 2 — v 2 scattering and final state corrections in 
order to accommodate the third unobserved particle radiated from the initial state 

S = (ZaPa+Pb) 2 = (P1+P2?, (4.78) 

* = (ZaPa-pif, (4.79) 

u = (z aPa -p 2 ) 2 . (4.80) 

In the limit of soft (p 3 = 0) or collinear (p 3 \\ p a ) particle emission, they satisfy the relation s + t + u = 
sz" — > for massless partons. 

The calculation proceeds as follows: We insert three additional <5-functions with respect to t, z", 
and z a into the general expression 

dPS(3) =/(2^n^|#^ L a + pb -± p \ (4.81) 

giving 



dPS (3) fm \dTT^Pi S ^) X d( , ^ 



dtdz"dz a J y ' AA (2vr) 

8(t - (z aPa - Pl f)5 (z" - ™ ) 5 (za - 22) . (4.82) 

V PaPbJ V PaPbJ 

Next, we integrate over the 5(pf) and the space-like components of the d-dimensional 5-function to 
eliminate p 2 - In the resulting expression 

dPS (3) f d d - 1 pxd d - 1 p 3 



I (2n)^2E 1 2E 2 2E 3 S [ Ea + Eh ~ g^ 



dtdz"dz a J (2n) 2d ~ 3 2E 1 2E 2 2E : 

6(t - (z aPa - Pl f)5 (z» - ») 5 L - ™* ) , (4.83) 

V PaPbJ V PaPbJ 

we now decompose p\ into its energy and angular components in the center-of-mass system of partons 
1 and 2 and p 3 in the overall center-of-mass system 

/o\ d — 3 d — 4 

, dPS l = f--^-, 2 -^-E d 1 - 3 dE 1 sm d - s 9de S m d -Ud<P-^^E* d - 3 dEI 

dtdz"dz a J (2ir) 2d ~ 3 8E 2 p (*=£) T (&=2) 

S m d ~ 3 x*dx*dcf>* 3 d I E a + E h - £ E j j 5(t + 2z a E a E 1 (l - cos 6)) 

5(^-25^(1-00.^))*^-^). (4.84) 

Integrating over the remaining 5- functions and the trivial azimuthal angle <\>\ up to 2ir, we arrive at 
- ■ ■■ — -z (ut) (l-z) ^ (l-Za-z) z a r sin 0— -, (4.85) 



dkk"dz a 7 1287r 3 r 2 (l-e) v ' v ; v a / a ■ r^ 

where iV^ is the normalization factor given in eq. ( [4.35| ). Finally, we can factorize this three particle 
phase space into 

dPS (3) =dPS (2) dPS w , (4.86) 
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where dPS^ 2 -* is the usual phase space for the two observed jets 1 and 2 from section 3.1 and 

is the phase space of the unresolved subsystem of partons a and 3. The integration measure is 

dm = izV -'^ (^M'f *■■*« (4.88) 

z a \l-ZaJ N^ r^(l-e) 

and the function 

H(z") = (1 - z")- 1+2e ( 1 - -?— ) =\ + 0{z") (4.89) 

can be approximated by 1 as it leads only to negligible terms of 0(y). The integration of dPS' r ' over 
z" G [0, —u/s], z a G [X a , 1], and <f> G [0, ir] is restricted to the singular region of z", when partons p a 
and j>3 are recombined, such that < z" < min{— u/s,y} = yj. 

4.2.5 Photon Initial State Corrections for Direct Photons 



Singularities from the initial direct photon arise from diagrams of type b) in figure 15. The photon 
here splits up into a quark- antiquark pair, and one of the two becomes a part of the photon remnant. 
If the (anti-)quark is collinear to the incoming photon, a simple 1/e pole is produced. Quadratic 1/e 2 
poles corresponding to soft and collinear singularities do not exist, since there is no direct coupling 
to a gluon line. In this specific case, it is therefore not necessary to perform a partial fractioning 
decomposition. 

The matrix elements for the 2^3 processes ^q — > qgg and ^q — > qqq are computed from the 
generic diagrams in figure O as in the case of final state corrections in section 4.2.2. Those for the 
processes with incoming anti-quarks and gluons are obtained by simple crossing of the diagrams or 
Mandelstam variables according to table 01. We keep only terms that are singular in the variable 

z" = ^ (4.90) 

PaPb 

and express the ten invariants for 2^3 scattering through z" , the longitudinal momentum fraction 
z a transferred from the photon to the parton entering the hard subprocess, and the variables s, t, and 
u defined in section 4.2.4: 

PaPb = 7T~ ( 4 - 91 ) 

PaPl = ^— (4-92) 

2z a s 

s — u 

PaP2 = (4.93) 

2z a s 

PaP3 = ^~z" (4.94) 

s —u 

PbPl = 7, z a (4.95) 

2z a s 

s -t , t . 

PbP2 = z a (4.96) 

2z a s 

PbP3 = ^-(1-Za) (4.97) 

^z a 

P1P2 = ^~z a (4.98) 

^z a 

P1P3 = ^- — 0--z a ) (4-99) 

2z a s 

s — u 

P2P3 = « (1-^a) (4.100) 

zz a s 
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Under these approximations, a singular kernel can be f actor ized out, which is universal for all processes 
of the type 76 — ► 123 and describes the splitting of a boson (in this case the photon) into two fermions 
(i.e. the quark-antiquark pair). Consequently, a parton from the photon scatters now off the parton 
in the proton and the relevant parton-parton Born matrix elements show up. In the case of 



\M\ 



2 



(s,t,u) 



q ^ sz" 



z 2 a + (l 



T, 



gg^gg 



(s,t,u) 



it is the quark-antiquark annihilation process into two gluons. For 



\M\ 
it is the process qq 

\M\ 



2 

■yq^qqq 



2 

ig^ggg 



s,t,u) 

. The result for 
(s,t,u) 



H sz" 



e 2 ey^ 



sz" 



£ + 0- 



4 + (l 



<■ qq->qq 



(s,t,u) 



T, 



qg-*qg 



(s,t,u) 



(4.101) 



(4.102) 



(4.103) 



is easily obtained by crossing (s <-> t) in the first process and multiplying by (—1). 

As the matrix elements do not depend on the variable <fi, the integration over <p in dPS^'-* is as 
trivial as before. In addition, the integration over the simple pole in z" can also be carried out easily. 
However, the integration over the momentum fraction z a involves a convolution with the photon 
spectrum in the electron. Since these exist only in a parametrized form too complicated to integrate, 
the ^-integration has to be done numerically and is written down explicitly in the final result 



dPS (r) |A4| 



7^123 ^) t ) 



s,t,u) 



dz a 



-e 2 9 2 ^ 



Pis 

2^ 



VV r (!- e ) 



X a 



r(l - 2e) 



'76^123 



[s,t,u). 



(4.104) 



The functions I^b- 


-,123 are given by 










ig^ggg\^i > ^) 


= 


\ - 1 P 
e2N c ? " 


-r(*«) + 2N c Pq * 


-j(z a )]n 


( l-*a N 


>+!] 


V z a , 


-l-yq^qqq{S, t, U) 


= 


\ - l P 
e2Nc q - 


-l(Za) + 2Nc P V 


_ 7 (z a )ln 


( 1-V 


i+f 


V Z a J 


ig^ggg\^i *■> ^v 


= 


\ - l P 
e2N c q ^ 


-lit*) + 2Nc Pq * 


_ 7 (z a )ln 


( l-z a ^ 


>+f 


K Z a J 



T, 



qq -+gg(s,t,U), (4.105) 

qq^qq{s,t,u), (4.106) 

T qg ^ qg (s,t,u). (4.107) 



T, 



Terms of 0{e) and 0{y) have been omitted as before. 

The collinear 1/e poles are proportional to the splitting function 

P q ^{z a ) = 2N c e 2 q P q ^ g {z a ) 



(4.108) 



for photons into quarks, where e q is the fractional charge of the quark coupling to the photon and 



P 



g^g 



'Za) 



1 



4 + (1 - Za)' 



(4.109) 



is the Altarelli-Parisi splitting function for gluons into quarks. This function appears in the evolution 
equation of the photon structure function as an inhomogeneous or so-called point-like term (see section 
2.3). Therefore, the photon initial state singularities can be absorbed into the photon structure 
function. The necessary steps are well known |53|, |68| . We define the renormalized distribution 
function of a parton a in the electron F a / e (X a , M%) as 



F a/e {X a ,Ml) 



1 dZa 

X a Z a 



a 



£ a7 <5(l - Z a ) + — R q ^ 1 (z a ,M a ] 

ZTT 



F, 



7/e 



X„ 



(4.110) 
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where R has the general form 

R a ^(z a ,M 2 a ) 



-P, 



<?<-7 



r(i- £ ) U^\ 

KZa) T(l-2e) I Ml ) +Uq ^ 



>o) 



(4.111) 



with C = in the MS scheme. The renormalized partonic cross section for 76 — > jets is then calculated 
from 

d<r(7& — ► jets) = da^b — ► jets) — — / d2: a i?„^_ 7 (2! a , M a )dcr(a6 — ► jets). (4.112) 

d<r(7& — > jets) is the higher order cross section before the subtraction procedure, and da(ab —> jets) 
contains the LO parton-parton scattering matrix elements T a b^i2(s,t,u). The factor Att^/M^ in 



eq. (4.111) is combined with the factor Att/j, / s in eq. (4.104) and leads to M* dependent terms of the 



form 



q^j{Z a ) 



4-ITfJ 2 



47T/X 2 

~mT 



<?^7 



(^a)ln 



Ml 



(4.113) 



so that the cross section after subtraction in eq. ( 4.112| ) will depend on the factorization scale M%. 



4.2.6 Proton Initial State Corrections for Direct Photons 

Initial state singularities cannot only show up on the direct photon side, but also on the proton side. 
The parton from the proton, which undergoes the hard scattering, will then radiate a soft or collinear 
secondary parton. A quark can, e.g., radiate a gluon as in figure p~5|c) , which will then not be observed 
but contributes to the proton remnant. As we can now have soft gluons, we find not only single but 



also quadratic poles in e. After singling out the matrix elements for the diagrams in figure 17 singular 
in the variable 

III PbP3 



Z 



PaPb 



(4.114) 



we therefore have to decompose them with partial fractioning using REDUCE [64]. The invariants 
are approximated quite similarly to the photon initial state corrections and are given by 



PaPb = 


S 

2z~b 


PaPl = 


s -t 
2zb s 


PaP2 = 


s —u 
2z b s 


PaP3 = 


2z b 


PbPl = 


s —u 
2z h s 


PbP2 = 


s -t 
2zb s 


PbP3 = 


S z'" 
2z h 


PlP2 = 


s 

2z h 


PlP3 = 


s -u 
- 2^ b - {l - Zb) 


P2P3 = 


s _t n ^ 

= 2^- {l - Zb) 



(4.115) 


(4.116) 


(4.117) 


(4.118) 


(4.119) 


(4.120) 


(4.121) 


(4.122) 


(4.123) 


(4.124) 
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The result for the first subprocess "yq — ► qgg turns out to be 

i 



\M\ 2 q ^ qgg (s,t,u) 



^V 8 -.,,, 



sz" 



4C F ((l-z b )(l-e)-2 + 



-2t/s 



z"> -t/s(l- z h ) 



-AN C 



-t/s 



-u/s 



z"> - t/s(l - z b ) z'" - u/s{\ - z b ) 



-yq^gq 



[s,t,u). (4.125) 



Since the initial quark has to couple to the photon, it cannot vanish in the beam pipe. The singular 
kernel only describes the radiation of a gluon from the quark, factorizing the leading-order QCD 
Compton process jq — > gq. According to table ||, we find the same result for incoming anti-quarks. 
In the four-fermion process 



\M\ 



2 _, 
iq^qqq 



2„2 4 ,6s 



j. \ l A 4 ' 

s,t,u) =e e q g [i 



1 



sz" 



1 + (1 - Zb) 

Zb 



£Z b 



ig^qq 



[s,t,u), 



(4.126) 



the initial quark necessarily produces a collinear final state quark. Thus, we have a non-diagonal 
transition of the quark-initiated process 73 — ► qqq into the gluon-initiated process 75 — ► qq. The other 
non-diagonal transition shows up in 



„2 2 4 ,6s 



\^\jg-*gqq,i («, t,u) = e'e*g*n' ,„ 



sz" 



C F 



(z 2 b + (1 - z b ) 2 )(l + e) - e 



T- 



■yq~+gq 



[s,t,u), 



(4.127) 



where the gluon-initiated process jg — » gqq is transformed into the quark-initiated process jq — ► gq. 
The splitting of a gluon into a collinear quark-antiquark pair is analogous to the splitting of the direct 
photon in the last section. However, the gluon also possesses a non-abelian coupling to other gluons, 
which can become soft or collinear in 



\M\ 



79^9<7<?.2 



[s,t,u) 



e 2 e^V fe 



+ 



sz' 1 
t/s 



-,N C 



2 1 — + z b {l-z h ) 
z b 



-u/s 



Z'" - t/s(l - Zb) z'" - u/s(l - z b ) 



r ig __gg(s,t,u) (4.128) 



and factorizes the photon-gluon fusion process 73 — ► qq. 

The complete list of proton initial state corrections given above has to be integrated over the 
phase space region, where parton 3 is an unobserved part of the proton remnant. The phase space 
for three-particle final states is taken from section 4.2.4, where (z" <-> z'") and (z a <-► z b ) have to be 
interchanged. The relevant integrals are calculated in appendix B. The result is 



dPS^\M\ 2 



76^123 



(s,t,u) 



dzb 2 2 2 Ae a s /4vr^ 2 \ r(l-e) 



A',, 



Zb 



ee « g » 2^ 



r(l - 2e) 



J, 



76^123 



(s,t,u), 



(4.129) 



where the functions J 7 6^i23 are given by 



.1 



iq-^qgg 



(s,t,u) 



2C F 



1 1 



Pq<-q(Zb) + 6(1 ~ Z bl l , , o _ 



1 1 / ~t\ 1 n -t n 

+ — 3 - 2 In — + - In 2 — + vr 2 
e* 2e V s J 2 s 



+l-z b + (l- z b ) In (yj^T^) + 2 ^+ (^ 
-2^^1nfl + ^i^ 



2& / 
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-2 In 



£ /l -z;, 

2fe 



1. « 1 

e t + 2 
I +21n 
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-2^1n(l+^i-^|+2-*-ln(l+^ii-^ 



1 - z b 



yjs z b 



l- Zb 



yjs z b 



J. 



-yq^qqq 



(s,t,u) 



"f<1^9<l\ S i i U )i 



--^rW**) + ^r p 9<-*^) ( ln ( yj hr ] + ' 



,1-^6 



•^6 



L ig-^qq 



(s,t,u), 



(4.130) 



(4.131) 



19~ > <W/>1 V^J '"> ^] 
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(4.132) 
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Here, we have used the function 
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In 1 + 



-u 1 — z b 



yjs z b 
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-ig^qq 



(s,t,u). (4.133) 



V 



I- z b 



(4.134) 
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as an abbreviation. As the integration over z b in eq. ( 4.129| ) runs from X b to 1, the +-distributions 
p9| ] are defined as 

D + [g}= [ dz b D(z b )g(z b ) - f dz b D(z b )g(l), (4.135) 



where 



x,, 



g(z b ) = —F b ,/ (— ) h(z b ), 
z b /y V z b 



(4.136) 



and h(z b ) is a regular function of z b [70]. This leads to additional terms not given here explicitly when 
eq. ( |4.135| ) is transformed so that both integrals are calculated in the range [X b , 1]. 

Some of the J 7 b^i23 contain infrared singularities oc l/s 2 , which must cancel against the corre- 
sponding singularities in the virtual contributions. The singular parts are decomposed in such a way 
that the Altarelli-Parisi kernels in four dimensions 



P q^q(Zb) 

Pg^q(zb) 
Pg^g(zb) 

Pq^g(Zb) 



C F 
C F 

2N C 
1 



1 + zl 3 r/ 



1 + (1 - z b ) 



Zb 
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zl + (1 - z b f 



+ — + z b (l -z b )-2 
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-N c - -N f 
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(4.137) 

(4.138) 

5(1 -Zb), (4.139) 

(4.140) 
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proportional to 1/e are split off. They also appear in the evolution equations for the parton distribution 
functions in the proton in section 2.2. The singular terms proportional to these kernels are absorbed 
as usual into the scale dependent structure functions 



F b/p (X b , M 2 ) = f 1 ** 
Jx h z b 



Zb 



"&<-&' 0*6) -^6 ) 



J- s R' 



Fb'/ P 



Xb 

Zb 



(4.141) 



5 bb >5(l 

where F b t/ p (X b /z b ) is the LO structure function before absorption of the collinear singularities and 

r(l-e) Avr/i 2 ^ 



R' b ^ b ,(z b ,Mi 



l -p, 

e 



'(zb) 



r(l - 2e) 



Ml 



+ C b ^ b ,{z b ) 



(4.142) 



with C = in the MS scheme. Then, the renormalized higher order hard scattering cross section 
d<r(7& — >jets) is calculated from 

da( 7 6 -» jets) = da(jb -» jets) - ^ f dz b R' b ^ b ,(z b , M^da^b' -» jets). (4.143) 

da^b — > jets) is the higher order cross section before the subtraction procedure, and da^b' — > jets) 
contains the lowest order matrix elements T jg - t gg(s,t,u) and T^ g ^ g g(s,t,u) in d dimensions. This 
well known factorization prescription [68, 71, ^(J removes finally all remaining collinear singularities. 
It is universal and leads for all processes to the same definition of structure functions if the choice 
concerning the regular function C" in ( 4.142| ) is kept fixed. Similar to the case of photon initial state 
singularities, the higher order cross sections in fl4,143| ) will depend on the factorization scale M b due 
to terms of the form P b ^ b i{z b ) ln(M 2 /s). 



4.2.7 Photon Initial State Corrections for Resolved Photons 

We now turn back to the case of resolved photons and consider their initial state corrections. We start 
with the singularities on the photonic side of the hard scattering cross section. For direct photons, 
there was only one possible singularity coming from the splitting of the photon into a quark-antiquark 
pair as shown in figure |l5|b). However, resolved photons contribute to the hard scattering like hadrons 
through their partonic structure. Therefore they produce similar poles to the proton initial state 
singularities as in figure |l~5|c) or those in section 4.2.6. These poles can be of quadratic strength due 
to the radiation of soft and collinear gluons off the quarks in the photon. 

We use the same O(a^) matrix elements that were already calculated for the final state corrections 
of resolved photons in section 4.2.3. The relevant diagrams are those for parton-parton scattering in 
figure [18|. As stated before, these diagrams show only the generic types and have also to be used 
in their crossed forms for the complete set of diagrams. All possible processes are listed in table |6] 
together with their matrix elements, initial spin and color averages, and statistical factors. Only the 
matrix elements singular in the variable z" = PaVz/PaPb are kept, where p% is the momentum of the 
parton soft or collinear to the original parton in the photon with momentum p a . Parton 3 will then 
be a part of the photon remnant. The matrix elements are expressed through the ten approximated 
invariants given in eqs. ( 4.91| )-( |4.100| ) and decomposed with the help of partial fractioning. It is then 
possible to factorize them into singular kernels and parts of the leading-order parton-parton matrix 
elements. 

Let us start with the process qq' - 
to the incoming quark a 

2 



qq'g in figure 18 a). Here, the final gluon can be soft or collinear 
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> w (s,t,u), (4.144) 
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Process 


Matrix Element \M \ 2 


qq 1 -> qq'g 


l\M\l^ M ( S ,t,u) + \M\l^ ag , Qa (s,t,u)]/[4N-£] 


qq' -> qq'g 


l\M\ 2 q^qqiq,i(u,t,s) + \M\ 2 qq ,^ qq , q ^u,t,s)]/[4N 2 } 


qq -> q'q'g 


[\M\ 2 qq ^ qq , q ^s,u)]/[w 2 ] 


qg -> qq'q' 


l\M\ 2 qg ^ qqW (s,t,u) + \M\ 2 g ^ qq , q ,(u,t,s) + \M\ 2 qg ^ qqlq ,(t,s,u) 
-\M\ 2 qq ^ qq , q , 2 (t,s,u)}/[8(l-e)N 2 C F ] 


qq -► qqg 


[\M{ 2 q/ ^ qq , g>l (s,t,u) + \M\ 2 qq/ ^ qq , gjl (s,u,t) + \M\ 2 qq ^ qqg (s,t,u) 

+2\M\ 2 qq ^ qqlq ^s,t,u)]/[AN 2 ]/2\ 


qq -► qqg 


+\M\ 2 qq ,^ qq , qa (u,t,s)]/[W 2 } 


qg -> qqq 


[\M\ 2 lql {s,t,u) + \M\ 2 , ql (s,u,t) + \M\ 2 qg ^ qqq (s,t,u) 

+\Mf q , q ,(u,t,s) + \M\ 2 qlql (t,u,s) + \M\ 2 qg ^ qqq (u,t,s) 

+\M\l g _ qq , q ,(t,s,u) + \M\\ g ^ qq , q ,{u,s,t) + \M\ 2 g ^ qqq (t,s,u) 

-2\M\ 2 q ^ M (t,s,u)]/[8(l - e)N 2 C F ]/2l 


qq -> 999 


[\M\ 2 qq ^ qqq (s,t,u)]/[AN 2 ]/3\ 


qg -> qgg 


[-\M\ 2 qq ^ qqq (t,s,u)/3 + \M\ 2 aa (s,t,u) + \M\ 2 qq ^ qqq , 2 (s,t,u)]/[8(l - e)N 2 C F }/2\ 


qg -*• qgg 


[~\M\ 2 qq _> qqq (t, u, s)/3 + \M\i^ Qq ^( U , t, s) + \M\ 2 qq _> qqq , 2 (u, t, s)]/[8(l - e)N 2 C F ]/21 


gg -> qqg 


[-\M\ 2 q ^ qqq At,s,u) + \M\ 2 qq ^ qqq (s,t,u)]/[W(l-e) 2 N 2 C 2 F ] 


gg ~> ggg 


l\M\ 2 qq ^ qqq (s,t,u)]/[l6(l-e) 2 N 2 C 2 F ]/3\ 



Table 6: Summary of 2 — ► 3 squared matrix elements for resolved photoproduction. 



so that the original quark will also participate in the Born process qq' — ► qq'. However, the quark 
can also go into the photon remnant. It will then radiate a gluon that scatters from the quark b with 
different flavor 
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(s,t,u) = g n 



1 
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For equal flavors in figure [l8|b), the interference contribution 
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(4.146) 



has the same kernel for the color factor (2C F — Nc), but a crossed version for the color factor C F , 
and is proportional to the Born interference contribution T qq ^ qq . If the final gluon is crossed into the 
initial state for the processes qg — > qq'q' and qg — > qqq with unlike and like quark flavors, it can split 
up into a collinear quark-antiquark pair quite similar to the splitting of direct photons in section 4.2.5. 
Consequently, the singular kernels in 
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L 99 ^99 
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'- qq^qq 



[s,t,u) 



(4.147) 
(4.148) 



are the same as in eqs. (|4.101|) - (|4.103 ), and we have the quark-quark scattering processes qq' — > qq' 
and qq — > qq on the tree level. The additional factor of (1 + e) is due to the averaging of initial gluon 
spins with 1/(2(1 — e)) ~ 1/2(1 + e), whereas the photons were averaged just by 1/2. The third process 
to be considered is taken from figure 18:). Obviously, only a gluon can go into the photon remnant 
so that in 
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we find the same simple pole for the color factor Cp as in eq. ( 4.144| ) and still have quark-antiquark 
annihilation into gluons in the Born process. The double poles are, however, more complicated and 
factorize only parts of the Born matrix elements. If a final gluon is crossed into the initial state, a 
soft or collinear gluon can also be radiated from the initial gluon 
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leaving a new kernel and a partly factorizable qg — > qg scattering process behind. Alternatively, the 
quark can go into the photon remnant 
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which leads to a gg — > gg Born process. For collinear quarks, only a single divergence is possible. The 
next kernel for gg — > qqg is already known from qg — > qq'q' and qg — > qqq 
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where a gluon splits into a quark-antiquark pair, but now with a different leading-order matrix element 
Q9 ~^ Q9- Finally in 
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it is clear that the Born process must also be completely gluonic and we can only have the gluon 
splitting into two gluons as in qg — > qgg. All these contributions have to be considered several times 
and also in crossed forms according to table ^. 

What remains to be done is the integration over the phase space of particle 3 in a region, where 
it can be considered a part of the photon remnant. We do so with the help of the phase space 
dPS' r ) calculated in section 4.2.4 and leave the ^-integration for numerics. Again, the reason is the 
analytically not integrable form of the photonic parton densities that have to be convoluted with the 
matrix elements. The result is 
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The integrated matrix elements are put into the functions I a b- 
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The absorption of the collinear poles 1/e proportional to the different Altarelli-Parisi splitting functions 
is handled in a completely analogous way as in section 4.2.6 for proton initial state corrections. The 
only difference is that the poles are absorbed into the photon structure function and not the proton 
structure function. As always, we omit terms of higher order in e and the invariant mass cut-off y. 

4.2.8 Proton Initial State Corrections for Resolved Photons 

Finally, the next-to-leading order O(a^) resolved photoproduction cross section receives also initial 
state corrections on the proton side. It is, however, not necessary to calculate these contributions 
again. The relevant diagrams in figure 18 are the same as in the last section as is the singularity 
structure in figure [l5|c). This is due to the fact that resolved photons behave like hadrons. 

The proton initial state formulas are obtained from those in section 4.2.7 by interchanging (z" <-> 
z'") and (z a <-> zt,), so that we consider matrix elements that are singular in the variable 



The parton b in the proton gives a fraction 
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(4.166) 



of its momentum to the 2^2 hard scattering process and the rest of (1 — z^) to particle 3 in the 
proton remnant. The list of approximated invariants is the same as in section 4.2.6 for proton initial 
state corrections of direct photons, and the list of contributing matrix elements is the same as in table 



The integration over the singular region of phase space gives 
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where the ^-integration is done numerically to allow for inclusion of the proton structure function. 
The functions 

Jab^l23(s,t,u) = I a b^l23(s,t,u) (4.168) 

are identical to those from the last chapter. For the IR singularities and finite contributions propor- 
tional to the 5(1 — Zft)-function, the integration can of course be carried out trivially. Thus, these 
singularities cancel against those from the virtual corrections for resolved photoproduction. The 
collinear singularities proportional to the Altarelli-Parisi splitting functions are now absorbed into the 
proton structure functions and not the photon structure functions as before. 
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4.2.9 Real Corrections for Direct 77 Scattering 

Real corrections to direct 77 scattering arise through the radiation of a gluon off one of the quark 



lines in the underlying Born process 77 — ► qq. This can be inferred from figure 19. The calculation 



Figure 19: 2^3 Feynman diagram for direct 77 scattering. 

of the final and initial state singular parts proceeds along the same lines as for direct and resolved 
photoproduction. 

Let us first consider the final state singularity. There, the approximated matrix element 
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is integrated over the singular phase space to give 



dPS^\M\ 2 ^ 123 (s,t,u) = e^^- 



i j ± £ a s ( 47T/u 2 \ T(l - e) 



2vr \ s 



r(l - 2e) 



■77-+123 



[s,t,u) 



with 



(4.170) 



3 2tH 



+ 7 - 21n 2 y F - 3lny F T, 



77^gg( s ) t, U). 



(4.171) 



Turning to the initial state, we find the approximated matrix element 
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The pole is proportional to the Altarelli-Parisi splitting function and is absorbed into the photon 
parton density. 



4.3 Finite Next-To-Leading Order Cross Sections 

We conclude this section with a summary of all singularities that appeared in the next-to-leading order 
cross section of jet photoproduction. There were three types of singularities: 

• UV singularities in the virtual corrections 
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• IR singularities in the virtual and real corrections 

• Collinear singularities in the initial state real corrections 

All of them were regularized dimensionally by going from four to d = 4 — 2e dimensions, where the 
regulator e had a positive sign for the ultraviolet (UV) and a negative sign for the infrared (IR) 
divergencies. 

The ultraviolet divergencies were encountered in the calculation of the virtual diagrams in section 
4.1. The diagrams had an additional inner "virtual" particle and were classified into self-energy 
diagrams, propagator corrections, box diagrams, and vertex corrections. The inner loop momenta 
could not be observed and had to be integrated up to infinity. The resulting UV singularities could 
be removed by renormalizing the fields, couplings, gauge parameters, and masses in the Lagrangian 
through multiplicative renormalization constants Zi, which show up in perturbation theory as counter 
terms order by order in the strong coupling a s . The counter terms were not given explicitly, so that all 
matrix element formulae in section 4.1 are already UV-divergence free and all fields, couplings etc. have 
to be considered physical and renormalized. As an example, the counter term for the QCD Compton 
graph 75 — > gq has the form 

\M\^ gq , CT (s,t,u) = W^)^ 

G + in i?) - t Nc ) r ^™ (s ' *• u) (4i75) 

in the MS scheme and leads to a logarithmic dependence of the cross section on the renormalization 
scale ii. 

The second type of singularities, IR divergencies, were produced at the lower end of the loop 
integration in the virtual corrections and through soft or collinear real particle emission in the real 
corrections. They were written down explicitly throughout the last sections and have to cancel ac- 
cording to the Kinoshita-Lee-Nauenberg theorem |6l]| . We demonstrate this cancellation separately 
for direct and resolved photoproduction in tables [7| and || and for direct 77 scattering in table ||[ 



Process 


Color Factor 


NLO Correction 


Singular Parts of Matrix Elements 


7<? -> gq 


c F 


Virtual Corr. 
Final State 
Initial State 




-£-I(3-2J(i)) 
+^ + i(3-2Z(t)) 
+^ + ^(3-2Z(t)) 


Tyq— ^gq\S) ^3 u ) 
Tyq—>gq(s, t, ll) 


N c 


Virtual Corr. 
Final State 
Initial State 




-^-^{f-2l(s) + 2l(t)-2l(u)) 

+ '(*)" *(«))" 
+i( - 1(8) + l(t)- l(u)) 


Tjq—>gq(s, t, ll) 

Tyq— >gq( s i ti u ) 
T'yq^gq(s, t, U) 


N f 


Virtual Corr. 
Final State 





Table 7: Cancellation of IR singularities from virtual, final state, and initial state NLO corrections 
for the direct partonic subprocesses and different color factors. 

There is only one generic 2 — ► 2 diagram 75 — > gq in direct photoproduction as shown in figure 
0, from which the photon-gluon fusion process can be deduced with the help of crossing relations. 
Therefore we show only this process in table |7[ but for the three contributing color factors Cp, Nc, 
and Nf . The real corrections for the first two classes come from the process jq — > qgg and are equally 
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divided among final and initial state singularities. The last class, however, occurs only in the splitting 
of the final gluon into an additional quark-antiquark pair with Nf flavors. This is the reason why no 
initial state singularity is present here. 

For resolved photoproduction, we have the four generic processes in figure ^. They are presented 
in table [8| and divided further into color factor classes. All other processes can be obtained through 
crossing. The real corrections for quark-quark scattering with different and like flavors arise simply 
from the emission of an additional gluon and factorize the complete Born matrix element. Final and 
initial state corrections contribute at equal parts. For processes involving more gluons, the situation is 
more complex. In qq — > gg, an additional gluon leads to different color factors depending on whether 
it is radiated in the initial state (Cf) or the final state (Nc)- A final gluon can also split up into 
Nf flavors accounting for half of the real corrections in the class Nf. The other half comes from the 
process qg — > qgg, where the initial gluon splits up into a quark-antiquark pair with Nf flavors. So far, 
the factorization property of qq — ► gg holds. However, the logarithmic contributions to the emission 
of a third initial or final gluon in the color classes 1, 8N^Cf, and 8NcCf only factorize parts of the 
leading order cross section. In the completely gluonic process gg — > gg it does not matter where the 
third gluon is radiated (color class Nc)- Still, a final gluon can split up into Nf flavors or an initial 
gluon can have come from Nf different quarks. Finally, the logarithmic contributions proportional to 
32NqCf only factorize parts of the Born cross section, but are symmetric under cyclic permutations 
of the Mandelstam variables as a completely symmetric process must be. 

For direct 77 scattering, there is only one Born matrix element 77 — * qq and only one color class 
Cf- The virtual singularities are canceled by the final state singularities alone as shown in table [|. 

The third and last class of singularities are those in the initial state from collinear real particle 
emission. Although they are generally classified as infrared singularities as well, they are not included 
in the tables [|], ||, and |^ above. These single poles proportional to the Altarelli-Parisi splitting functions 
P 9 «_ 9 (z), Pg<-. q (z), Pqi-g(z), and Pg^ g i z \ do not cancel against similar poles from virtual corrections. 
They are absorbed into the renormalized photon and proton structure functions according to the MS 
scheme and leave behind a logarithmic dependence of the hard cross section on the factorization scales 
M a and M b . 

The finite next-to-leading order cross section for the photoproduction of two jets was given in 
section 2.5 as 

d^cr dcr 
dE 2 dmdm =Y^x a F a/e {x a ,Ml)x b F h / p {x b ,Ml) — {ab^p lP 2) (4.176) 

with the partonic cross section 

dcJ , , x 1 ^TPT dPs(2) /„ ^ 

^^ P1P2 ) = -\M\^. (4.177) 

We can now return to the physical four dimensions letting e — * everywhere in the calculation. The 
strong coupling constant a s and the parton densities in the electron F a / e (x a , M%) and in the proton 
F b / p (xb, M b ) are renormalized and defined in the MS scheme. 

5 Numerical Results for Photoproduction 

In this section we present numerical results first for single and dijet NLO cross sections in complete 
photoproduction, as they were defined in section 2.5. We use the analytical results that have been 
calculated in leading order in section 3 and in next-to-leading order in section 4. All UV, IR, and 
collinear initial state singularities canceled or could be removed through a renormalization procedure, 
leading to finite results in four dimensions. Equivalent numerical results for single and dijet NLO cross 
sections in complete 77 collision are presented in the next chapter. 

The cross section formulas are implemented in a FORTRAN computer program, which consists 



of four main parts as explained in table 10. The first part contains the analytical results of sections 
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Process 



Color Factor 



NLO Correction 



Singular Parts of Matrix Elements 
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4Z(«) - 2Z(t))' 



Virtual Corr. 
Final State 
Initial State 



+\(4l(s)-2l(u)-2l(t)) 
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4N C C F + %-\^- 16N C C F ^ 



Virtual Corr. 
Final State 
Initial State 



+ 7 P (l " 

-h m t 



Virtual Corr. 
Final State 
Initial State 



1 (u 




+i(i + i)m+Ku)) 
+iii±iim±M) 
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7V C 



Virtual Corr. 
Final State 
Initial State 



_ 4 _ 22 

3s 

, 2 , 11 
+ ^ + 37 



'gg^gg 



1 99^99 



L gg-+gg 



(s,t,u) 
(s,t,u) 
(s,t,u) 



Virtual Corr. 
Final State 
Initial State 



+ -L Tgg^gg ( s > u ) 

2 T~1 



q,l—gg\S, I . u) 

(s,t,u) 



S2N C C F 



Virtual Corr. 
Final State 
Initial State 




^gg^gg 



2H + 



r,tU I t +U 

Z 7Z + "Wr 



+ cycl. perm.^ 
+ cycl. perm. 



(Z(s) 3 - 2$ + ^ + cycl. perm. 



Table 8: Cancellation of IR singularities from virtual, final state, and initial state NLO corrections 
for the resolved partonic subprocesses and different color factors. 
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Process 


Color Factor 


NLO Correction 


Singular Parts of Matrix Elements 


77 -> 99 


c F 


Virtual Corr. 
Final State 




_£_l( 3 -2J(t)) 
+£ + J(3-2J(t))" 


T~/y-+qq(s , t, 1l) 



Table 9: Cancellation of IR singularities from virtual and final state NLO corrections for direct 77 
scattering. 



Part No. 


Number of Jets 


Contributions 


1 


2 


Analytical contributions in LO and NLO 


2 


2 


Numerical contributions in jet cone 1 


3 


2 


Numerical contributions in jet cone 2 


4 


3 


Numerical contributions outside jet cones 



Table 10: Organization of the FORTRAN program into four main parts. 



3 and 4, and the numerical parts are needed for an implementation of the Snowmass jet definition 
(see section 5.1). Obviously, the third part is only needed for two-jet cross sections. For one-jet cross 
sections, this region of phase space is included in part four. 

The main task of the computer program is to integrate over different parameters. For example, 
total or partly differential cross sections are integrated over Ex, rji and/or r/2- Furthermore, the 
momentum fractions of the initial photon in the electron x a , the parton in the photon y a , and the 
parton in the proton Xj, have to be integrated. Additional momentum fractions z a and Zb appear in 
the initial state corrections. Finally, the numerical contributions have to be integrated over the phase 
space of the third parton Er 3 , and 03. The two 5- functions for the momentum fractions of the 
partons going into the hard 2^2 scattering X a and X\» 

' (*. - w e E E r.*-*) - d s ( x » - w P £ ") ■ <"> 

reduce the total number of integrations by two. The non-trivial task of computing up to seven- 
dimensional integrals is solved with the Monte Carlo routine VEGAS written by CP. Lepage |ff^] , 
which adapts the spacing of the integration bins to the size of the integrand in the bin. 

The input parameters for our predictions will be kept constant throughout this chapter, if not 
stated otherwise. All cross sections are for HERA conditions, where electrons of energy E e = 26.7 GeV 
collide with protons of energy E p = 820 GeV. Positive rapidities rj correspond to the proton direction. 
We use the Weizsacker- Williams approximation of eq. ( |2.13| ) with Q^ax = 4 GeV 2 as in the ZEUS 
experiment, but do not restrict the range of longitudinal photon momentum x a in the electron. For 
the parton densities in the proton, we choose the next-to-leading order parametrization CTEQ3 in the 
MS scheme, which already includes HERA deep inelastic scattering data. The corresponding A value 
of A^ 4 ) = 239 MeV is also used in the two-loop calculation of the strong coupling constant a s with four 
flavors. We do not use the one-loop approximation for leading order calculations, because the effects 
of the next-to- leading order hard scattering contributions are then isolated. The parton densities in 
the photon are taken from the NLO fit of GRV and transformed from the DIS 7 into the MS scheme. 
The renormalization and factorization scales are equal to Ex- We use the Snowmass jet definition with 
R = 1, no jet double counting, and no R scp parameter. This makes the cross sections independent 
of the phase space slicing parameter y, which is fixed at y = 10~ 3 . This value is sufficiently small to 
justify the omission of 0{y) terms in the analytical calculation (see section 5.1). 
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In section 5.1, we check our analytical results with respect to y-cut independence and with two 
different existing programs for direct and resolved one-jet production. Section 5.2 contains studies 
of the dependence of the cross sections on renormalization and factorization scales and on various 
cancellation mechanisms. Theoretical predictions for one- and two-jet cross sections are presented in 
sections 5.3 and 5.4, before we conclude this section with a comparison of our calculation to data from 
the HI and ZEUS collaborations in section 5.5. 



5.1 Check of the Analytical Results 

The first check of our analytical results for the photoproduction of jets in next-to-leading order of QCD 
has already been given in section 4.3. There it has been shown that all infrared divergencies cancel 
in a consistent way between the virtual and the real corrections, giving a finite cross section in four 
dimensions. We have therefore missed none of the singular terms. This cancellation mechanism could 
only work since we integrated the 2^3 matrix elements over regions, where two final state particles 
or an initial and a final state particle had an invariant mass Sjj = (pi + Pj) 2 smaller than a fraction 
Vf,i,j of the center-of-mass energy s. Naturally, this leads to a dependence of the cross sections on 
these phase space slicing parameters uf,i,j- We choose y cu t = Uf = Hi = Vj in the following. 

One possibility to deal with the dependence of the cross section on the invariant mass cut-off 
y cut is to use it as a definition for the experimentally observed jets. As we noted already in section 
2.4, the experimental jet definition has only a correspondence in theory beyond the leading order. 
Furthermore, we mentioned a special kind of algorithm, i.e. the JADE cluster algorithm, which uses 
exactly the invariant mass criterion to cluster hadrons into jets. For e + e~-experiments, one can then 
identify the theoretical (partonic) with the experimental (hadronic) cut-off 

theory experiment (K 0\ 

Vcut — Vcut ■ 

In this case, the corresponding cross sections are for exclusive dijet production. 

The numerical value of y cut is constrained in a two-fold manner. First, the real corrections contain 
simple and quadratic logarithmic terms In y cut , In 2 y cut , coming from the single and quadratic 1/e 
poles, which force the cross section to become negative for y cut < 10~ 2 . Therefore one should choose 
Vcut > 10~ 2 . Second, we have always assumed the singular region to be small in the real corrections. 
This means we have omitted terms like y C utlny C ut, Vcut, ••• This forces us to take y cut <C 1. The 
typical value for the JADE algorithm is therefore given by y cu t — 10~ 2 . 

The omitted terms of 0(y C ut) can be calculated numerically to improve the precision of the cal- 
culation. This makes it necessary to include the full and un-approximated 2^3 matrix elements in 
the FORTRAN program, which are then integrated numerically in the region between the analytical 
cut-off 2/cut° ry an d the experimental cut-off y^{ 3crimcnt . It is important to note that the full 2^3 
matrix elements have to be integrated after partial fractioning. Regions, where a single invariant is 
small but is not a pole in the matrix element, cannot be neglected. The theoretical invariant mass cut 
is then reduced to a purely technical variable, on which the physical prediction will no more depend. 
If one integrates over the phase space outside the experimental cut-off as well, one arrives at inclusive 
dijet cross sections. In this case, the full 2^3 matrix elements contribute as a leading-order process 
through the production of three jets, where the third jet is unobserved. 

At HERA, we have to deal with a partly hadronic process, for which cluster algorithms are not so 
well suited. We will therefore follow the two experiments HI and ZEUS and use the Snowmass cone 
algorithm, where hadrons i are combined into a single jet, if they lie inside a cone in azimuth-rapidity 
space from the jet center. The experimental cone size R was already defined in section 2.4 



Ri = V Xm - vj) 2 + - <Pj) 2 < R (5.3) 

and will be chosen to be R = 1 in the following. As R takes the role of ^ experiment ; we expect the sum 
of the analytical two-body and numerical three-body cross sections to be independent of y^ ^ ■ This 
is the second decisive test of our analytical results. 
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Figure 20: Inclusive single-jet cross section (f a / dEx dq for direct photons at Et = 20 Ge V and n = 1, 
as a function o/y cu t- The analytical (dashed) and numerical (dotted) contributions have a quadratic 
logarithmic dependence on y C ut, which cancels in the full next-to-leading order result (dot-dashed curve). 
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Figure 21: Inclusive single-jet cross section cPa/dExdr] for resolved photons at Et = 20 GeV and 
7] = 1, as a function of y cut . Again the sum of two-body (dashed) and three-body (dotted) curves is 
y C ut-independent (dot-dashed curve) like the leading order (full) curve. 
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In figure ^ we plot the inclusive single-jet cross section d 2 a/dET<lr] for direct photons at Et = 
20 GeV and ij = 1 as a function of y cu t- The leading order prediction (full curve) is trivially independent 
°f Vcut- The analytical two-body contributions (dashed curve) exhibit a quadratic logarithmic depen- 
dence on y cut and turn negative below y cut = 2.3 • 1CT 2 . We add the numerical contributions inside 
and outside the jet cone into the three-body contributions (dotted curve), which are also quadratically 
logarithmically dependent on y C ut, but positive. The sum of two- and three-body contributions is the 
full inclusive next-to- leading order result (dot-dashed curve) and is independent of y cut . This proves 
that the y cu t-dependent finite terms in our analytical results are correct. Towards very small values of 
2/cut — 10~ 5 , the dot-dashed curve drops slightly. This is due to the limited accuracy in the numerical 
integration of the three-body contributions and can be remedied with increased computer power. 

Figure 21 is the analogue to figure 20 for resolved photoproduction. The two-body curve turns 
negative at a slightly smaller value of y cut = 1.5 • 10 -2 . As the number of contributing partonic 
subprocesses is much larger than in the direct case, more computer time was needed for similar 
statistical accuracy. 

The ultimate test for our analytical and numerical results consists in the comparison of our pre- 
diction with existing predictions, that were obtained by Bodeker []53| (direct case) and by Salesch 
[ff^ l (resolved case) with different methods. This can only be done for the single-jet inclusive cross 
section dPa/dExdrj, since dijet cross sections are not available. Bodeker's results are obtained with 
the subtraction method of Ellis et al. [58] to cancel soft and collinear singularities. Salesch calculated 
finite cone corrections in addition to the results by Aversa et al. [74], who used the so-called small 
cone approximation. Here, a small jet cone radius 5 is used to slice the phase space and isolate the 
final state divergencies. The dependence on 5 cancels, when the finite cone corrections are added. We 
take exactly the same parameters in our calculation and in the reference calculations. Furthermore, 
we choose y cu t = 10 -3 to ensure that we are not sensitive to the omission of terms of 0(y cu t). 

The Et dependence is checked in figure 22 with Bodeker for direct photons and in figure ^3| with 
Salesch for resolved photons at r/ = 1. Our prediction (full curve) agrees with the older calculations 
(open circles) at a level better than 1%. The numerical statistics are so good that the error bars in 
Bodeker's and Salesch's programs are not seen. For comparison, we include the leading order dotted 
curves, where the predictions are identical. 

We compare the same inclusive single-jet cross section as a rapidity distribution for Et = 20 GeV. 
We choose now a linear scale for the ordinate in figures 24 for direct and 25 for resolved photopro- 
duction. This makes it possible to see the error bars in the reference calculations (open circles) which 
are at the same level as the agreement between the different programs. Errors of ~ 1% present in our 
program are not shown. With more computer time, the agreement would certainly be even better. 

From these comparisons, we conclude that also the cut-independent finite terms in our analytical 
results and in our numerical FORTRAN program are correct. Even more, it was possible to check the 
transverse energy, rapidity, and other distributions not shown here for every single direct and resolved 
partonic subprocess. We found perfect agreement everywhere. 

It is clear that the results presented in this work can also be applied to the calculation of inclusive 
jet cross sections in pp collisions by just replacing the photon parton distributions with those of the 
antiproton and considering only the resolved case. We could have extended our tests further by 
comparing with published results for the pp case. Such results have been given in 5£] using the 
subtraction method and in |)(| based on the phase-space slicing method. In connection with work on 
the large transverse momentum production of single inclusive jets in pp collisions pTj, in which the 



program based on [73] was used, we compared with results from [^] and found good agreement. A 
comparison with the two-jet results in |)6| has not been attempted. 

After completion of this work, two new complete calculations of jet production in low Q 2 ep 
collisions have been presented. In p8j| , the subtraction method for canceling infrared and collinear 
singularities was applied to obtain various inclusive one- and two-jet cross sections. These authors 
compared inclusive single-jet and two-jet cross sections with our published results [12| and found good 
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Figure 22: Inclusive single-jet cross section (Pa/ 1 dExdrj for direct photons at r\ = 1, as a function of 
Et- The subtraction method results (open circles) agree with our phase space slicing prediction (full 
curve). 
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Figure 23: Inclusive single-jet cross section cPa/dETdn for resolved photons atrj = l, as a function of 
Et- The small cone approximation and finite cone corrections (open circles) agree with our invariant 
mass cut method (full curve). The leading order (dotted) curve is shown for comparison. 
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Figure 24: Inclusive single- jet cross section (fa/dErdr] for direct photons at Et = 20 GeV, as a 
function of r\. Our predictions (full curve) agree with Bodeker's predictions (open circles). 
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Figure 25: Inclusive single-jet cross section Sa/dETdn for resolved photons at Et = 20 GeV, as a 
function of n. Good agreement is seen for Salesch's (open circles) and our calculation (full curve). 
The leading order (dotted) curve is shown for comparison. 
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agreement [98|. Another work [99] uses also the phase-space slicing method, but with two separate 



parameters for the infrared and the collinear singular regions as in their earlier work [52]. Their results 



compare to our results in [12|. A direct thorough comparison of the results obtained with the same 
input has not been done yet. 

5.2 Renormalization- and Factorization-Scale Dependences 

From the numerical checks in section 5.1, we can now be sure to have a reliable computer program for 
the calculation of jet cross sections in photoproduction. We will use this program here to study the 
dependence of these cross sections on the three relevant scales: 

• the renormalization scale \i 

• the factorization scale for the photon M a 

• the factorization scale for the proton Mj, 

We expect a reduced dependence in next-to-leading order for all three scales with respect to the leading 
order. As stated in sections 4.1, 2.2, and 2.3, we take //, M a , and Mf, to be of O(Et)- Therefore, 
we will always plot the single-jet inclusive cross section d 2 a/dETdrj as a function of scale/ Et in the 
following. Direct and resolved photoproduction will be presented separately at a fixed transverse jet 
energy of Et = 20 GeV and at a rapidity of rj = 1, where the cross sections are at their maximum. 

We start with the dependence on the renormalization scale fi. Leading order 0(aa s ) cross sections 
depend on fi only through the running of the strong coupling constant a s (/i 2 ), which has to be 
implemented in the one-loop approximation 

as tf) = ™L (5.4) 

(33-2iV»ln£ 



for consistency. The results are shown as dotted curves in figures 26 and 27. For direct photons, the 
curve drops by one half when going from \x = 1/4Et to jj, = AEt- The resolved photon hard cross 
section is one order higher in a s (0(a 2 )) and therefore drops even by a factor of four. This strong 
l/ln()U 2 )-behavior is slightly weakened in the two-loop approximation 



aV) = 1* 6(163- 19^) h-Q.^! 

(33 -21V,) In £ I (33-2JV»» ]njj> 1 



where one adds further logarithmic terms in \i , These curves are shown in dashed form in figures ^ 
and 27. However, going from one-loop to two-loop a s affects not so much the shape but the overall 
normalization. The cross sections are reduced by 20% in the direct and by almost 40% in the resolved 
case. Although not physically relevant, the dashed curves can be compared to the full curves and 
enable us to see solely the effect of the explicit logarithmic terms 

iffHdVii+talHd, (5.6) 

e V s J e s 

that appear in the one- loop corrections to the hard cross section in section 4.1. They cancel partly the 
a s -dependence, so that the full curves in figures ^ and ^ are much flatter compared to the leading 
order curves with one- and two- loop a s . The direct curve varies by only 20% and exhibits a maximum 



near \i)Et — 1/2. According to the principle of minimal sensitivity [75], fi = \/2Et would then be 
the optimal scale. The dependence of the resolved photoproduction cross section is also weakened 
and amounts to a factor of two in NLO. It meets the one-loop LO result at fi ~ \/2Et- This scale is 
sometimes preferred since the NLO corrections are small in its vicinity. 
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Figure 26: Inclusive single-jet cross section (fa / ' dExdq for direct photons at Ex = 20 GeV and n = 1, 
as a function of /i/Et- From fi = 1/AEt to fi = AEt, the leading order predictions drop by one half. 
The two-loop curve (dashed) is 20% smaller than the one-loop curve (dotted). The next-to-leading 
order (full) curve only varies by 20%. 
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Figure 27: Inclusive single-jet cross section cPa/dExdr] for resolved photons at Et = 20 GeV and 
■n = 1, as a function of /jl/Et- From \i = 1/AEt to = 4Et, the leading order predictions drop by 
a factor of four. The two-loop curve (dashed) is 40% smaller than the one-loop curve (dotted). The 
next-to-leading order (full) curve only varies by a factor of two. 
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The next scale that we will study is the factorization scale in the photon M a . It separates the soft 
parton content in the photon given by the parton distribution function F a /^(y a , M 2 ) from the hard 
partonic cross section da a b(ab — ► jets). As discussed in section 2.3, direct and resolved processes are 
only separable in leading order. The direct parton distribution is then simply a 5- function 5(1 — y a ) 



and does not depend on any scale. This can be seen from the dotted curve in figure |28|. The resolved 
part depends, however, strongly on the scale M a as could already be inferred from e + e~-data on the 

We expect the M a -dependence to be dominated by the asymptotic pointlike 



-Fr/ structure function 
contribution to 



a 



1 



-a(x B ) + b(x B ) 



(5.7) 



«s(M2) 

that behaves like In M 2 . Clearly we see this behavior in the leading order dot-dashed curve in figure 
which is identical to the dotted curve in figure ^ and is obtained with the GRV photon distributions. 
As the GRV distributions are fitted in the DIS 7 -scheme, we transform them into the MS scheme 
by adding the correct finite terms independent of M a . The renormalization scheme dependence was 
studied by Bodeker, Kramer, and Salesch |76j and will not be studied here again. The main result is 
that the dependence disappears in next-to-leading order, because the transformation of the resolved 
part is compensated through scheme dependent terms in the initial state singularities of the direct 
photon. The renormalization of the singularities in section 4.2.5 



F ah {y ai M 2 a ) 
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through the 1/e pole 
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4-KfJ 2 
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M 2 



(5.10) 



We can see this negative logarithmic behavior (—hiM 2 ) in the dashed curve in figure 28 for next- 
to-leading order direct photoproduction. It has the opposite sign than the LO resolved contribution, 
so that the sum of both will eventually be almost independent of M a (full curve in figure |28| ). The 
next-to-leading order resolved contribution has a similar, but slightly steeper shape than the leading 
order and is shown as a full curve in figure |2^. Therefore, the dependence of the complete NLO 



calculation on the photon factorization scale is slightly stronger than in the full line of figure 28. The 
difference will be compensated in next-to-next-to-leading order (NNLO) of direct photoproduction, 
which is beyond the scope of our calculation. These findings agree with the earlier analysis of the 
M a dependence of the inclusive single-jet cross sections in |76| and constitute a nice check of our 
formalism. 

For the proton, we have a similar scale Mb as in the photon case that separates the soft and the 
hard part of the hadronic and partonic cross sections. Of course, the proton is not point-like and does 
not have a direct component. The scale dependence of the proton distribution functions F b j p (x} 1 ,M 2 ) 
is governed by the AP equations. If we solve them iteratively, we can assume that in 



dF b/p (x b ,M 2 
dlnM fe 2 



(o) 
as 



dz b 



P, 



b^b' 



1 ^b/p 



(Zb) 
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Figure 28: Inclusive single-jet cross section d 2 a / dE? drj for direct photons at Et = 20 Ge V and n = 1, 
as a function of M a / Et- The LO direct curve (dotted) is independent of M a as is the sum (full curve) 
of NLO direct (dashed) and LO resolved (dot-dashed) contributions. 
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Figure 29: Inclusive single-jet cross section cPa/dETdr] for resolved photons at Et = 20 GeV and 
7] = 1, as a function of M a / Et- The logarithmic behavior of the resolved part is slightly steeper in 
NLO (full curve) than in LO (dotted curve). 
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the starting distributions Ffty(z b ) and the coupling constant are scale independent. Then, the 
integral over d In M? can be carried out leading to a simple logarithmic behavior 



F b/p (x b ,M^) 



a 



(o) 



i 



2ir J z b 

x b 



x b \ „(0) 

7 b ) V 



(5.12) 



For the LO resolved photon contribution, this approximation is obviously good enough as can be seen 
from the dotted curve in figure 31. The cross section drops by 20% when going from M b = 1/AEt 
to M b = 4Et- Of course, the full solution of the AP equations is not so simple, and the LO direct 



photon contribution in figure 3C deviates slightly from a simply logarithmic graph. The variation of 
the cross section here amounts only to 6.5%. In next-to- leading order, initial state singularities arise 
as in the photon case and are absorbed into the proton structure function 



F b/p (x b ,M^) 
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Again, finite terms accompany the pole in the MS scheme 
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which depend on M b through 



-Pb^b'(Zb) 



ATTfJ 2 



(5.15) 



These logarithms cancel the leading logarithmic behavior of the parton distribution functions in the 
proton. Consequently, the full next-to-leading order curves in figures ^ and [H] depend only very 
weakly on the renormalization scale M b in the proton. Finally, we consider the total scale dependence 



for complete photoproduction in figure 32. The dependence of the inclusive single-jet cross section on 
M = fx = M a = Mb is dominated by the renormalization scale dependence in the one-loop (dotted 
curve) and two-loop (dashed curve) approximation for a s . Since resolved photoproduction is more 
important than direct photoproduction at Et = 20 GeV, figure ^ resembles figure more than 
figure The total scale dependence is reduced by almost a factor of two in next-to-leading order 
(full curve). 



5.3 One-Jet Cross Sections 

In the last section, we have studied the dependence of the inclusive single-jet cross section on the 
various scales, which cannot be predicted from theory. It was shown that the next-to-leading order 
cross section is much less dependent on the renormalization and factorization scales than the leading 
order cross section, due to cancellation mechanisms. We continue to study one-jet cross sections here, 
but shift the emphasis towards experimentally relevant distributions. The two main experimental 
observables in the double differential cross section 

are the hardness or transverse energy Et of the observed jet and its orientation in the detector, i.e. the 
angle 9 it forms with the beam axis. Instead of this angle, we use the pseudorapidity rj = — ln[tan(0/2)] 
as usual. As the electron is only weakly deflected, it stays in the beam pipe like the proton remnant. 
Therefore, the jets are homogeneously distributed in the azimuthal angle (f>, and we have integrated 
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Figure 30: Inclusive single-jet cross section (fa/dExdr] for direct photons at Et = 20 GeV and r\ = 1, 
as a function of Mb/Ex- The weak behavior of the LO (dotted) curve is even more reduced in NLO 
(full curve). 



over <j) in the theoretical prediction. The single-jet cross sections for direct photoproduction have 
already been published ||, those for resolved photoproduction are shown here for the first time with 
our phase space slicing method. Using a different method similar results have been published in 



73, 78, 7 



We start with distributions in the jet transverse energy Et, which is at the same time a measure 
for the resolution of partons within the photon and the proton. In figure 33, we plot the direct photon 
contribution for transverse energies between 5 and 70 GeV. In this interval, the cross section drops 
by almost six orders of magnitude. The next-to-leading order curves are always above the leading 
order predictions. To separate the curves for three different rapidities, we multiplied those for r] = 
and r) = 2 by factors of 0.1 and 0.5, respectively. At r] = 0, the phase space restricts the accessible 
range in Et to Et < 50 GeV. The so-called fc-factor, i.e. the ratio of NLO to LO, drops from 1.8 at 
5 GeV to 1.2 at 70 GeV for r/ = 1, so that the higher order corrections are more important at small 
scales, where the strong coupling a s is large. Had we calculated the leading order with one-loop a s , 
the correction would, however, not be so large. 

The corresponding curves for resolved photoproduction are shown in figure 34. Here, the curves 
for r] = 1 and r/ = 2 lie closer to each other in spite of the rescaling for r/ = 2. This already hints at a 
broader maximum in the rapidity distribution compared to the direct case (see below). The NLO and 
LO curves lie further apart due to the large /c-factor of about 1.8 over the whole £/r-range. This can 
be understood from the non-perturbative nature of resolved photons, where higher order corrections 
are more important at all scales due to the many channels involved. 

The rapidity distribution for direct photons is presented in figure |3^ at a transverse energy of 
Et = 20 GeV. In the available phase space of rj € [—1,4], the cross section exhibits a rather broad 
maximum and steep edges. The /c-factor is 1.25 in the central region. 

We have already mentioned that the corresponding distribution for resolved photons is expected 
to have an even broader maximum. This proves to be true in figure [if], where the maximum is also 
slightly shifted in the proton direction from n = 1 to rj = 1.5. The reason can be found in the lower 
longitudinal momentum of photon components with respect to a direct photon. The next-to-leading 
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Figure 31: Inclusive single-jet cross section d 2 a / dEx dq for resolved photons at Et = 20 GeV and 
11 = 1, as a function of M^/Et- The approximately logarithmic dependence in LO (dotted curve) is 
reduced in NLO (full curve). 



order cross section is about 80% larger than the leading order cross section, which is compatible to 
the observation of a large /c-factor from the £/r-distribution above. 

We now look at the interplay of direct and resolved photoproduction in inclusive single-jet cross 
sections. Figures [37| and |3^ present the leading order and next-to- leading order predictions for complete 
photoproduction (full curves), which are the sums of the direct (dotted) and resolved (dashed) curves 
already presented above. The rapidity of the observed jet is fixed at rj = 1. Obviously, the resolved 
photon component is dominating at low transverse energies due to the photon structure function 
and unimportant at large transverse energies, where the perturbative point-like coupling is large. In 
leading order, the intersection of the resolved and direct curves lies near 26 GeV and is shifted towards 
37 GeV in next-to-leading order. The resolved process has larger NLO corrections than the direct 
process, so that the region, in which the latter is important, moves to higher Et- 

The rapidity distributions for full photoproduction are shown as full curves in figures H and [|j 
for leading and next-to-leading order at a transverse energy of Et = 20 GeV. One expects the direct 
photon to be important mostly in the electron direction at small or negative rapidities. This behavior 
can be seen in both figures, but only close to the boundary of phase space below rj = 0. 

In figures 41 and 42, we plot the dependence of the single-jet inclusive cross section on the jet 
cone size R in the Snowmass convention pq| . As discussed in section 2.4, a theoretical definition for 
jets containing more than one parton is only possible in next-to-leading order. Therefore only the 
full curves depend on R, whereas the leading order curves are constant. The dependence has the 
functional form of the next-to-leading order result [73, 77] 



, = a + b\nR + cR\ (5.17) 
dETdn 

For the direct photon contribution in figure f|l], the leading and higher order predictions are equal at 
R ~ 1 for one-loop a s (dotted) and at R ~ 0.6 for two-loop a s (dashed). For the resolved photon 
contribution in figure the situation is different. The stable points lie at R ~ 0.7 for one-loop 
(dotted) and R ~ 0.12 for two-loop a s . 
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Figure 32: Complete inclusive single-jet cross section (Fo-/dE T dri at E T = 20 GeV and rj = 1, as a 
function of M/Et for M = [A = M a = M&. The strong dependence in LO (dotted and dashed curves) 
is reduced in NLO (full curve). 



5.4 Two-Jet Cross Sections 

We now turn to two-jet cross sections, where one does not integrate over the second rapidity 772 or, 
alternatively, over the momentum fraction x a of the parton in the electron. The differential cross 
section 

d^cr dcr 
dE 2 dmdm =^2 x aF a / e (x a ,M^)x b F b/p (x b ,M^) — (ab^p 1 p 2 ) (5.18) 

then yields the maximum of information possible on the parton distributions and is better suited to 
constrain them than the observation of inclusive single-jets. Since dijet production is a more exclusive 
process than one-jet production, the cross sections are smaller and require higher luminosity or longer 
running time in the experiments. This is the reason why HI and ZEUS have only recently started 
to analyze dijet data and why we can present here the first theoretical calculation for complete dijet 
photoproduction. 

It is important to note that only in leading order the transverse energies of the two observed 
jets balance (E^ = Ej> 2 = Ex). In next-to-leading order inclusive cross sections, there may be a 
third unobserved jet which must have full freedom to become infinitely soft. Therefore, the transverse 
energies of both jets cannot be observed without spoiling infrared safety. This is an artifact of fixed 
order perturbation theory and will go away in O(aa^), where one may as well have a fourth jet. We 
will calculate similar distributions as in the last section, i.e. in the transverse energy of the first jet, 
Er 1; and in both observable rapidities rji and r\i. Yet, the rapidities r\\ and 772 always belong to the 
two jets with largest transverse energy in the event. As in the one-jet case, the direct distributions 
have already been presented in a recent paper ||, whereas the resolved and complete cross sections 
are shown here and partly in 



First, we look at the distributions in the transverse energy E^ in figure |43| for direct photons. We 
fix r/i at 771 = 1 and 772 at three different values of 772 = 0, 1, 2. The curves for 772 = and 772 = 2 are 
rescaled by factors of 0.1 and 0.5 as before. We expect the dijet cross sections to be smaller than the 



single-jet cross sections in figure |33|. Indeed, for Et, Et x = 5 GeV and 77, 771, 7/2 = 1, the cross section 
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Figure 33: Inclusive single-jet cross section Sa/dExdn for direct photons as a function of Et for 
various rapidities r) = 0, 1, 2 in LO and NLO. The cross section for rj = (rj = 2) is multiplied by a 
factor of 0.1 (0.5). 
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Figure 34: Inclusive single-jet cross section d 2 a/ ' dE^dn for resolved photons as a function of Et for 
various rapidities rj = 0, 1,2 in LO and NLO. The cross section for r] = (rj = 2) is multiplied by a 
factor of 0.1 (0.5). 
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Figure 35: Inclusive single-jet cross section (Pa/dErdr) for direct photons as a function of rj for 
E T = 20 GeV in LO and NLO. 
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Figure 36: Inclusive single-jet cross section d 2 o~ / dExdn for resolved photons as a function of rj for 
E T = 20 GeV in LO and NLO. 
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Figure 37: Inclusive single-jet cross section d 2 a / dE^dn for complete photoproduction at rj = 1, as 
a function of Et- The full curve is the sum of the LO direct (dotted) and LO resolved (dashed) 
contributions. 
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Figure 38: Inclusive single-jet cross section dFa/dExdri for complete photoproduction at n = 1, as 
a function of Et- The full curve is the sum of the NLO direct (dotted) and NLO resolved (dashed) 
contributions. 
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Figure 39: Inclusive single-jet cross section d 2 a/ dExdr] for complete photoproduction at Et = 20 Ge V, 
as a function of n. The full curve is the sum of the LO direct (dotted) and LO resolved (dashed) 
contributions. 
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Figure 40: Inclusive single-jet cross section d 2 a/ dErdw for complete photoproduction at Et = 20 Ge V, 
as a function of rj. The full curve is the sum of the NLO direct (dotted) and NLO resolved (dashed) 
contributions. 
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Figure 41: Inclusive single-jet cross section d 2 a/ dE^dn for direct photons at Et = 20 GeV and rj = 1, 
as a function of the jet cone size R. Only the NLO (full) curve and not the LO curves with one- 
(dotted) or two-loop (dashed) a s depends on R. 




Figure 42: Inclusive single-jet cross section d 2 a/ dExdn for resolved photons at Et = 20 GeV and 
rj = 1, as a function of the jet cone size R. Only the NLO (full) curve and not the LO curves with 
one- (dotted) or two-loop (dashed) a s depends on R. 
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drops by almost an order of magnitude from 10.6 nb to 1.29 nb. The /c-factors in dijet production 
are basically the same as in single-jet production. There is, however, one novel feature: for 772 = 0, 
the third jet in next-to-leading order opens up some additional phase space, so that Er 1 can go up to 
Et x < 55 GeV instead of only 39 GeV in leading order. 



A similar behavior is seen for resolved photoproduction in figure 44. At the lowest value of 



Et 1 = 5 GeV and back-to-back jets (rji = t]2 = 1) 5 the two-jet cross section (4.91 nb) is even more 
than one order of magnitude smaller than the one-jet cross section (67.5 nb). The ratio of NLO to 
LO is 1.7 for resolved photons over the whole E?± -range. The phase space for 772 = increases from 
Et x = 39 GeV to Er x = 50 GeV. This slightly smaller value is due to the reduced center-of-mass 
energy in resolved photoproduction, as some of the energy always goes into the photon remnant. 

Next, we present the dependence of the cross sections on the two rapidities in form of the three- 
dimensional lego-plots ^5] and 46. The leading order is always shown on the left side and is completely 
symmetric in 771 and 772. For the next-to- leading order cross sections on the right hand sides of figures 
[45| and 4€ , this is no longer exactly true due to the presence of a "trigger" jet with transverse energy 
£/Ti> which is fixed at E^ x = 20 GeV. The next-to-leading order lego-plots are only approximately 
symmetric. This can best be seen at the bottom of the contour plots, where at least one of the two 
observed jets is far off the central region. The NLO predictions are considerably larger than the LO 
predictions, especially in the resolved case. 

This becomes even clearer when we plot the projections of the lego-plots for fixed rji = 0, 1, 2 
and 3. In figure ^?], we plot the leading and next-to-leading order distributions in 7/2 for direct 
photoproduction. It is clearly seen that the second jet tends to be back-to-back with the first jet, 
since the maximum always occurs at 772 — 771. However, at 771 = 3 this is no more permitted by phase 
space. We obtain the same A:-factor of 1.25 in the central regions as in single-jet production. 

The 772-distributions for resolved photons in figure 48 are considerably broader than those in figure 
[47] due to the smearing of the hard cross sections with the distribution function of partons in the 
photon. The maxima of the plots are also not so much dominated by kinematics but more by the 
quark and gluon structure of the photon in different x regimes. They do not lie at 772 = 771 any more. 
Therefore, dijet rapidity distributions are best suited to constrain the photon structure. We will come 
back to this point in the next section 5.5, when we compare similar plots to data from ZEUS. The 
/c-factors range from 1.65 in the central regions to more than 3 in the proton forward direction. The 
shapes of the distributions are very similar in LO and in NLO. The absolute values make, however, 
an important difference. 

Direct and resolved contributions are now added to give physical, complete photoproduction re- 
sults. The dijet cross section is first plotted as a function of the transverse energy E^- Figure f49| 
gives the LO result, figure 50 the NLO result. Like the direct and resolved cross sections alone, the 
full two-jet cross sections are about an order of magnitude smaller than the one-jet cross sections (see 
figures 37 and 38). The point where direct and resolved contributions are equally important is lowered 
towards Er 1 = 20 GeV in leading order and Et x = 30 GeV in next-to-leading order, so that direct 
photons are better observed in dijet production. 

If one plots the complete two-jet cross sections as a function of 772, the different behaviors of direct 
and resolved photons add up to the full curves in figures 51 and [5^. These plots are best suited to 
decide in which rapidity regions one can look best for the resolved photon structure. We have already 
seen that this will be in situations where the two jets are not back-to-back, e.g. for 771 = and positive 



772 in the upper left plots of figures 51 and 52. On the other hand, the proton structure can best 
be studied with direct photons, when the cross section is not folded with another distribution. A 
possible scenario is 771 = and negative values of 772. This is especially interesting for the small- x 
components of the proton like the gluons and the quark sea. Another interesting observation is that 
the relative importance of direct and resolved processes changes dramatically when calculating dijet 
photoproduction in next-to- leading order O(aa^): resolved processes are much more important at 
Etx = 20 GeV than one would have guessed from a leading order estimate. 
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Figure 43: Inclusive dijet cross section S"oj dET x dr\\dr\2 for direct photons as a function of Et x for 
rji = 1 and three values of 7/2 = 0, 1,2. The cross section for 7/2 = (rj2 = 2) is multiplied by 0.1 (0.5). 
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Figure 44: Inclusive dijet cross section cfa/ dE^drjid^ for resolved photons as a function of Et x for 
rji = 1 and three values of 7/2 = 0, 1,2. The cross section for 772 = (7/2 = 2) is multiplied by 0.1 (0.5). 
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Figure 45: Inclusive dijet cross section d?a/ dEr^dqidrfz at Et x = 20 GeV for direct photons, as a 
function 0/771 and 772 . The LO plot (left) is exactly symmetric, the NLO plot (right) only approximately. 




Figure 46: Inclusive dijet cross section cfia/ dE^drjid^ at Et x = 20 GeV for resolved photons, as a 
function ofr\\ andr]2- The LO plot (left) is exactly symmetric, the NLO plot (right) only approximately. 
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Figure 47: Projections of the LO (full curves) and NLO (dashed curves) triple differential dijet cross 
section for direct photons at Et x =20 GeV and fixed values 0/7/1 = 0, 1, 2, and 3, as a function of 




Figure 48: Projections of the LO (full curves) and NLO (dashed curves) triple differential dijet cross 
section for resolved photons at Et x =20 GeV and fixed values of 7/1 = 0, 1, 2, and 3, as a function 
0/7/2. 
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Figure 49: Inclusive dijet cross section £ 'a / ' dE^dijid^ for full photoproduction at n\ = r]2 = 1 as 
a function of Et x - The full curve is the sum of the LO direct (dotted) and LO resolved (dashed) 
contributions. 
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Figure 50: Inclusive dijet cross section d?a / dE^dnid^ for full photoproduction at r/i = 772 = 1 as 
a function of Et x - The full curve is the sum of the NLO direct (dotted) and NLO resolved (dashed) 
contributions. 
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Figure 51: Projections of the complete triple differential dijet cross section at = 20 GeV and fixed 
values 0/771 = 0, 1, 2, and 3, as a function 0/772. The full curve is the sum of the LO direct (dotted) 
and LO resolved (dashed) contributions. 




Figure 52: Projections of the complete triple differential dijet cross section at Et x = 20 GeV and fixed 
values 0/771 = 0, 1, 2, and 2>, as a function of '772. The full curve is the sum of the NLO direct (dotted) 
and NLO resolved (dashed) contributions. 
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5.5 Comparison of Photoproduction Results to HI and ZEUS Data 

In this section we compare the next-to-leading order calculation to recent one- and two-jet data from 
the HI and ZEUS collaborations at HERA. Both collaborations have continuously measured various 
cross sections for the photoproduction of jets since HERA started running in 1992. With the increased 
luminosity in recent years, the data have improved and many aspects of jet production could be studied. 
In our comparison we restrict ourselves to the measurements of one- and two-jet cross sections just 
recently published which are based either on 1994 or 1995 data. In particular, we shall compare with 



the inclusive single-jet data of 1994 from the ZEUS collaboration [82], with the inclusive dijet data of 



1994 from ZEUS [S3], with inclusive dijet data of 1995 from ZEUS [83], and with inclusive dijet data 
of 1994 from the HI collaboration @. 

We start with the single-jet cross section d 2 a/drjdET integrated over Et > -Er min as a function of r\. 
This cross section da/drj has been measured in the r/ range between -1 and 2 and with the Et thresholds 
ET min = 14, 17,21, and 25 GeV. The cross section da/drj for Et > 14 GeV has also been measured 
in three different regions of W: 134 GeV < W < 190 GeV, 190 GeV < W < 233 GeV, and 233 GeV 
< W < 277 GeV. The measurements refer to jets at the hadron level and are performed for two cone 
radii in the rj — (ft plane, R = 0.7 and R = 1 using the iterative cone algorithm PUCELL. The complete 
data have already been compared to our next-to-leading order calculations in the ZEUS publication 
|B2]. Therefore, we show only a selection for some specific kinematical ranges and compare them with 



the data. Our results for da/drj for Et > 17 GeV, R = 1 and R = 0.7 are shown in figures 53 and 



respectively, and are compared to the ZEUS data. The error bars in figures |53j, |5j, |55|, and |56| only 
contain the statistical error. The systematic error and the uncertainty associated with the absolute 
energy scale is not included, which adds an additional 30% error (see |82]). The theoretical predictions 



include resolved and direct processes in NLO. For the proton, the CTEQ4M 1 31 ] parton densities have 



been used. For the photon distribution, the GRV-HO ||33|j, converted to MS factorization, and as an 



alternative set the recent GS96 [38] parametrizations have been chosen as input. The renormalization 



and factorization scales have been put equal to Et, and a s was calculated with the two-loop formula 
with A^jg = 296 MeV as used in the proton parton densities. In figures [5^ and [54|, two curves are 
presented for both photon distribution sets labeled as R sep = 2R and R sep = 1R. They correspond 
to two choices of the -R sep parameter. Since our calculations include only up to three partons in the 
final state, the maximum number of partons in a single jet is two. Therefore, the overlapping and 
merging effects of the experimental jet algorithm cannot be simulated in the theoretical calculation 
[pq ]. To account for these effects, the i? se p parameter was introduced [85]. It has the effect that two 



partons are not merged into a single jet if their separation in the rj — (j) plane is more than R scp . Then 
i? scp = 2R means that no further restriction is introduced and the cone algorithm is applied in its 
original form. Experimentally, the two extreme values of R sep = 2R and 1R correspond to a fixed cone 
algorithm (like EUCELL) and to the Ut clustering algorithm (like KTCLUS), whereas an iterative 
cone algorithm (like PUCELL) is described by some intermediate value. In both calculation and data 
analysis, the maximal virtuality of the photon is equal to Qj^ax = 4 GeV 2 , and the full W range, which 
corresponds to 0.20 < x a < 0.85 in the EPA formula, is used. Looking at figures |5^ and [54|, we observe 
that the behavior of the measured cross sections is different for R = 0.7 and R = 1. For R = 1, the 
shape of the cross section is well described for -1 < rj < 0.5. For higher values of rj, the data stay 
almost constant as a function of rj, whereas the theoretical curves decrease as a function of rj for both 
radii R = 1 and R = 0.7. However, when R = 0.7 is used, the shape and magnitude of the NLO results 
agree quite well with the measured differential cross section in the entire ry range. This is also the case 
in the comparison for the lower Et threshold, -Er min = 14 GeV, with R = 0.7 shown in ||82|| . For the 
higher Et thresholds, 21 and 25 GeV, the NLO predictions give a good description of the measured 
cross sections in magnitude and shape for both cone radii R = 0.7 and R = 1 (see |8^] ). In general, the 
choice R = 0.7 should be preferred for the comparison between data and theory. Figure shows that 
the predictions with the GS96 parametrization of the photon parton distributions agree better with 
the data than the GRV-HO parametrization which is above the data for both R se p parameter values 
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Figure 53: rj dependence of the inclusive single-jet photoproduction cross section integrated over Et > 
17GeV with jet cone size R = 1. We compare our NLO prediction with GRV and GS96 photon parton 
densities and the two extreme R sep values to 1994 data from ZEUS. 
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Figure 54: r] dependence of the inclusive single-jet photoproduction cross section integrated over Et > 
17GeV with jet cone size R = 0.7. We compare our NLO prediction with GRV and GS96 photon 
parton densities and the two extreme R sep values to 1994 data from ZEUS. 
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over the whole rj range. Concerning the R sep parameter the curve for i? S cp = 1-R is in somewhat better 
agreement than the curve for i? se p = 2fl. We have checked that a value of -R SC p = 1.4R gives the best 
agreement. This supports a recent study of the jet shape function, which depends sensitively on this 
parameter [B6|. By comparing with recent measurements of this jet shape by the ZEUS collaboration 
|]8^1 , it was found that i? sep = 1.5R gives very good agreement with the jet shape data for PUCELL 
in the same r\ and Et range |86| . 

A comparison of da/drj for Et > 14 GeV in different W regions has also been presented. As 
an example, we show da/drj as a function of rj for the largest W range: 233 GeV < W < 277 GeV 
(corresponding to 0.55 < x a < 0.85) in figure |55| (R = 1) and in figure |56] (R = 0.7). Whereas the 
R = 1 theoretical cross section (figure |55|) agrees for low values of 77, it disagrees in the high 77 region 



with the data. This disagreement shows up particularly in the high W range [82]. The measured 
differential cross section is again well described by the NLO calculation for R = 0.7 (see figure ^).Q 

The excess of the measured cross section with respect to the calculations in the high r\ range and 
for the smaller Et thresholds for R = 1 is supposed to be due to the presence of the underlying 
event in the data which yields a larger amount of extra energy lying in the jet cone for R = 1 but 
not for the smaller cone R = 0.7. From figure 53 it is clear that the excess occurs only in the large 
r\ range where the resolved cross section dominates and where additional interactions of the photon 
and proton remnants are supposed to occur which are not included in the NLO calculations. These 
deviations between NLO theory and the data at large r\ and smaller Et's were found earlier [KM] when 
the theoretical predictions were compared with the 1993 ZEUS data ||. 

In conclusion, we can say that the NLO calculations describe reasonably well the experimental 
inclusive single-jet cross section for jets defined with R = 0.7 in the entire r\ range and for R = 1, if 
Et is large enough. 

Next, we compare the NLO predictions with inclusive dijet cross sections measured by ZEUS [p3|] 
and HI f84|| . Inclusive two-jet cross sections depend on one more variable as compared to the inclusive 
one-jet cross sections considered above. Therefore they are supposed to give a much more stringent 
test of the theoretical predictions than the inclusive one-jet cross sections. For the comparison with 
the data it is essential that in the theoretical calculations the same jet definitions are introduced as in 
the experimental analysis. Furthermore it is important that the theoretical calculations contain the 
same cuts on the kinematical variables as for the measured cross sections. First experimental data for 
inclusive two-jet production have been published by the ZEUS collaboration in || and |89|]. The more 



recent ZEUS analysis based on the 1994 data taking presented in [80] and recently in [^J] extends the 
earlier analysis in ^ based on 1993 data in several ways. The larger luminosity obtained in 1994 lead to 
a reduction of the statistical errors as well as allowing for the measurement of the cross section at higher 
Et, a region, where uncertainties due to hadronization of partons into jets and of underlying event 
effects are reduced making the comparison with the NLO predictions more meaningful. Furthermore, 
the ZEUS collaboration applied three different jet definitions: two variations of the cone algorithm 
@ called "EUCELL" and "PUCELL" , and the A; T -cluster algorithm "KTCLUS" as introduced for 
hadron-hadron collisions |9Cj. The two cone algorithms treat seed finding and jet merging in different 
ways. Since the NLO calculations contain only up to three partons in the final state, these experimental 
seed finding and jet merging conditions cannot be fully reproduced. This ambiguity is largely reduced 
in the /cr-cluster algorithm. In the NLO calculations, the two cone algorithms can be simulated by 
introducing R SC p already considered for the one-jet cross sections. The EUCELL definition corresponds 
to R sep = 2R, whereas PUCELL is best simulated with -R scp = 1.4R for the Et range considered in the 
experimental analysis (see above) |8(], |4q ]. The /c^-duster algorithm for hadron-hadron collisions is 
identical to using -R se p = R- Therefore by introducing the i? S ep parameter into the NLO calculations 
of the two-jet cross sections, all three jet finding definitions used in the experimental analysis can 
be accounted for. The ZEUS results of the earlier analysis []8^] have been compared to the NLO 

1 The wiggles in the curves in figures |55| and are due to insufficient numerical accuracy and have no physical 
significance. 
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Figure 55: r\ dependence of the inclusive single-jet photoproduction cross section integrated over 233 
GeV < W < 277 GeV with jet cone size R = 1. We compare our NLO prediction with GRV and 
GS96 photon parton densities and the two extreme R sep values to 1994 data from ZEUS. 




Figure 56: r] dependence of the inclusive single-jet photoproduction cross section integrated over 233 
GeV < W < 277 GeV with jet cone size R = 0.7. We compare our NLO prediction with GRV and 
GS96 photon parton densities and the two extreme i? sep values to 1994 data from ZEUS. 
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predictions in |T^] for the case of the KTCLUS algorithm and in |l0| for the EUCELL algorithm. 
These comparisons were done for the differential cross section da/df}, where r\ = 1/2(771 + 772) is the 
average rapidity of the observed jets with Et larger than ET min for both observed jets. This common 
cut on the Et of both jets causes some theoretical problems as has been noticed already some time 
ago @. 

The new measurements [83] are for the triple differential cross section d? a / dExdriid^ using the 
/cr-cluster or the PUCELL algorithm. The jet with the highest Ej- (leading jet, Et = Et x ) is required 
to have Et > 14 GeV and the second highest-lir jet to have Et 2 > 11 GeV. This cross section 
is symmetrized with respect to 771 and 772 and therefore double counted. By this symmetrization, 
the experimental ambiguity of determining the leading jet is avoided and the measured cross section 
corresponds to the calculated cross section where Et corresponds to the trigger jet which is not 
necessarily the leading jet with the highest Et- In order to have a handle to enhance direct over 
resolved photoproduction, one determines also the variable || [9l]] 



„obs _ Ziffle" 
2x a E e 



x obs = (5.19) 



where the sum runs over the two jets of highest Et and x a E e is the initial photon energy. x^ BS 
measures the fraction of the photon energy that goes into the production of the two hardest jets. The 
LO direct and resolved processes populate different regions of x® BS : x® BS = 1 for the direct process 
and x® BS < 1 for the resolved process. In NLO, the direct process populates also the region x® BS < 1. 
To obtain a measurement of the enriched direct photoproduction cross section, the cut x^ BS > 0.75 
is usually introduced. 

Figure [57] shows da/dET for six independent regions in the (7/1,7/2) plane. The upper curves 
in each plot give the dijet cross section for the entire x® BS region, and the lower curves which, to 
separate the curves, all have been scaled down by a factor of 5, present the cross section for the 
direct 7 region x!^ BS > 0.75. The NLO calculations are performed for Q^ax = 4 GeV 2 and 134 GeV 
< W < 277 GeV using the same parton densities for the proton and the photon as in the inclusive 
one-jet calculations. Furthermore we use R se p = 1.4i? with R = 1 to simulate the PUCELL algorithm. 
We compare with the corresponding data from ZEUS p3| analyzed with the PUCELL algorithm. 
The agreement between the data and the theoretical predictions is quite reasonable. Except near 
the backward regions (the last two 771,772 regions), the cross sections for the GRV and GS96 photon 
densities are very similar. To discriminate between them, the experimental errors must be reduced. 
We observe that the cross section for x® BS > 0.75 increases as compared to the full cross section (all 
x° BS ) as Et increases in agreement with the prediction from the calculation. This is the effect of the 
direct component which shows in general a flatter distribution with increasing Et than the resolved 
cross section [ |l~2"| , |10[ . We emphasize that the magnitude as well as the shape of the measured cross 
section is well reproduced by the calculations except for the first Et bin and for the case that both 
jets are in the region — 1 < 77^2 < 0. In this region (the last plot), the predictions lie above the data. 
The same cross section has been calculated also for the /cr-cluster algorithm. The results for the GS96 
photon densities have been presented together with the corresponding experimental data in [83]. The 



agreement between data and theory is quite similar. In |83y also data for da/dr]2 for Et > 14 GeV 
in three regions of 771 and the KTCLUS algorithm are compared to our NLO calculations with the 
GRV-HO and GS96 parametrizations of the photon densities. Both shape and magnitude of the cross 
sections are roughly reproduced by the calculations except for the small 771 region. Comparisons for 
the same cross section with the PUCELL cone algorithm have been done, too, but are not shown here 

& 

The comparison between measurements and calculations of the triple differential cross section 
shown so far covers the dependence of this cross section in all three variables Et, ??i and 772. Another 
equivalent set of variables consists of the dijet invariant mass Mjj, the rapidity of the dijet system yjj, 
and the scattering angle 9* in the dijet center-of-mass system. The dijet invariant mass is obtained 
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Figure 57: dependence of the symmetrized dijet photoproduction cross section integrated over 
different rapidity bins. We compare our NLO prediction with GRV and GS96 photon parton densities 
and the full and upper range of x!? BS to preliminary 1995 data from ZEUS. 
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from the relationship 

Mjj = 2E Tl E T2 [coshfai - m ) - cos(0! - 2 )] , (5.20) 

where and (f>2 are the azimuthal angles of the two jets in the HERA frame. For two jets back-to-back 
in <j) and with equal Ej<, 

Mjj = 2E T cosh [(771 - r/ 2 )/2] = 2E T /sm9* (5.21) 

and cos 9* = tanh[(r/i — 772) /2]. For events with more than two jets, the two highest £/r jets are used 
to calculate Mjj. 

The distribution in the dijet mass Mjj provides an additional test and is sensitive to the presence 
of resonances that decay into two jets. The dijet cross section as a function of cos 9* is sensitive to 
the parton-parton dynamics in the direct and resolved contributions |52|1 . Direct processes involve 
quark propagators in the t and u channels leading to a characteristic angular dependence proportional 
to (1 — I cos#*|) _1 . In the case of the resolved process, t-channel gluon exchange processes dominate 
which lead to an angular dependence proportional to (1 — | cos 9*\)~ 2 . This exchange rises more steeply 
with increasing | cos 9*\ than in the case of the direct processes. This different behavior in the angular 
dependence for resolved and direct processes was observed for Mjj > 23 GeV [p9|] . We have calculated 
the NLO cross section dcr/dcos #* for x° BS < 0.75 (= "resolved") and x° BS > 0.75 (= "direct") and 
have confirmed the different behavior in the two x^ BS bins as shown in the data (not shown here). 

The cross sections da / dMj j and da / dcos 9* have been measured recently using the sample of dijet 
events found with the PUCELL cone and the KTCLUS cluster algorithms. This analysis also includes 
the 1995 data which have an even higher statistics than used in the previous analysis (8t| based on 1994 
data. These cross sections have been measured in the kinematic region Q 2 < 4 GeV 2 , 0.2 < x a < 0.85 
as used previously. The two jets with highest Et are required to have Er x , Et 2 > 14 GeV and the 
rapidities of these two jets are restricted to —1 < 771,7/2 < 2.5. The cone radius is R = 1. The cross 
section da/dMjj has been measured in the Mjj range between 47 GeV and 120 GeV integrated over 
I cos 9*\ < 0.85. The cross section dcr/dcos 6>* has been measured in the interval < |cos6>*| < 0.8 
integrated over Mjj > 47 GeV. The experimental results @ ||] for da/dMjj are shown in figure EE 
separately for the two jet definitions, where the cross sections for the Izt algorithm have been scaled 
down by a factor of 10 in order to separate the data obtained for the two jet definitions. Systematic 



errors are available, but only statistical errors are shown here. In figure |58j, two NLO curves are 



compared to the measurements, the full curve being the result for the GRV-HO and the dashed curve 
the result for the GS96 parametrization of the photon densities. For the lower Mjj values, the GRV 
density seems to describe the data better than the GS96 density. However, we have to consider that 
the data points have an additional systematic error from the energy scale uncertainty |83||, which is 



also not shown in figure |58|. It is remarkable that the predictions for both jet algorithms, PUCELL 
and KTCLUS, agree well with the data over the full range of Mjj where the cross sections exhibits a 
fall-off of almost three orders of magnitude. 

The cross section dcr/dcos 8* as a function of | cos 6*\ between and 0.8 is plotted in figure 59. Here 
only the results for the PUCELL algorithm are shown. For the KTCLUS algorithm, the corresponding 
cross section is presented in [|g3[| . Again two curves are shown, for the GRV and GS96 photon densities, 
respectively. The theoretical curves agree reasonably well in magnitude and shape with the data. If it 
were not for the systematic error and the energy scale uncertainty, the comparison of data and theory 
would lead to a preference of the GRV over the GS96 density. We conclude, that the NLO calculations 
account reasonably well for the shape and the magnitude of the measured da/dMjj and dcr/dcos 9* 
as well as the triple differential cross section d 3 cr/dl!/rdr/id772 as a function of Et for various bins in 
7/1 and 772. 

The last comparison between data an our NLO theory concerns the double-differential inclusive 
dijet cross section d 2 cr/d2^ BS d-ET published just recently by the HI collaboration [B4|. This HI 
analysis is based on 1994 data. The photoproduction events have been selected with the constraint 
<2max = 4 GeV 2 and 0.2 < x a < 0.83. The jets were constructed with the cone algorithm with cone 
size R = 0.7. The implementation of the cone algorithm in the HI analysis is using a fixed cone and is 
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Figure 58: Mjj dependence of the dijet photoproduction cross section integrated over | cos 6**| < 0.85. 
We compare our NLO prediction with GRV and GS96 photon parton densities to 1995 data from ZEUS 
taken with the PUCELL and KTCLUS jet algorithms. 
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Figure 59: | cos#*| dependence of the dijet photoproduction cross section integrated over Mjj > 47GeV. 
We compare our NLO prediction with GRV and GS96 photon parton densities to 1995 data from ZEUS 
using the PUCELL jet algorithm. 
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therefore similar to the EUCELL algorithm used by ZEUS. The rapidities of all jets are restricted to 
the region —0.5 < rj < 2.5. In this specific two-jet analysis, Et is the average transverse energy of the 
two jets with the highest Et- The average rapidities of the two jets were between < (rji + r/2) /2 < 2 
their difference being \rji — 772 1 < 1, which corresponds to | cos 9*\ < 0.46. These two cuts on 771 and 772 
ensured that the jets are in a region with good measurements of the hadronic energy in the detector. 
The transverse energies of the jets were restricted further to the range 

'f-Tf 2 ' < 0.25, (5.22) 



and Et was required to lie above 10 GeV. This cut and the cut ( |5.22 ) ensure that the transverse 



energy of both observed jets is above 7.5 GeV (to avoid underlying event problems) without using the 
same Et cut for both jets, which would cause problems in the NLO calculations |8l]]. The observable 
x obs wag ca l cu l a ted from the same formula (|5.19| ) as in the ZEUS analysis. 

The NLO calculations of the inclusive dijet cross section d 2 cr/dx^ >BS dlog 10 (ii ; ^/GeV 2 ) is based 
on the same parton distributions for the proton and photon, respectively, as used in the previous 
sections. The scales are chosen equal to Et as in the comparisons above. We have chosen R scp = 2R, 
which we believe simulates best the fixed cone algorithm in the HI analysis. Otherwise the same cuts 
on Et 1 ,Et 2 ^'Hi^2 are applied as in the experimental analysis of the two-jet data. Our predictions 
are shown in figure again for the GRV-HO (full line) and the GS96 (dotted line) photon parton 
distributions. The curves are compared to the measured cross sections using the statistical and 
systematic errors added in quadrature. The data and the theoretical cross sections in the various 
x^ BS bins between 0.1 and 1 have been multiplied by factors 10™(n = 0, 1, 2, 6) in order to separate 
the results for the seven bins in x® . Our calculations give a good description of the data in magnitude 
and shape as a function of Et and x^* BS , except for some deficiencies in the two highest x® BS ranges. 
The GRV-HO and the GS96 parton distribution functions each give a satisfactory description of the 
measured cross sections with slight preference for GS96 in the lowest bin 0.1 < x® BS < 0.2. It is 
remarkable that even for x^ BS < 0.3 the theoretical cross section agrees so well with the data. This 
region was always a problem in earlier analyses due to the underlying event problems (see for example 
the inclusive single-jet cross sections for R = 1 of the ZEUS collaboration discussed above). These 
problems are apparently avoided in this analysis by choosing R = 0.7 and applying sufficiently large 
Et cuts. The deviations of the calculated from the measured cross sections could easily be due to 
hadronization effects which are not included in the NLO calculations. 

The cross section d 2 a/dx° BS dlog 10 (^/GeV 2 ) clS EL function of x® BS for fixed Et should show 
more clearly the dependence on the photon parton distribution sets. Therefore, we have calculated 



this cross section for Et = 11 GeV and 13 GeV. The result is shown in figures 61 and 62 for the GRV 



HO and the GS96 parametrizations and compared to the HI data [84], again including statistical and 
systematic errors added in quadrature. The agreement with the data is reasonable for both Et's. It 
seems that the hadronic maximum at x® BS ~ 0.3 is reduced from Et = 11 GeV to Et = 13 GeV, 
whereas the pointlike/direct peak at x® BS > 0.8 is enhanced flattening out the valley in between. 
This might be due to the QCD evolution of the VMD part of the photon structure function and 
the increased importance of the anomalous piece and the direct contribution at larger scales. For 
x° BS < 0.3, we see a somewhat better agreement for the prediction with GS96 which we already 



noticed in connection with the comparison in figure |60j. Both parton distribution sets give very similar 
predictions in the intermediate x^ BS range. The different results in the low and high x^* BS regions 
show us clearly how the experimental data must improve in accuracy before we can discriminate 
between the two alternative parametrizations of the photon distributions used in our calculations. 
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Figure 60: Er dependence of the symmetrized dijet photoproduction cross section integrated over 
different x® BS bins. We compare our NLO prediction with GRV and GS96 photon parton densities 
and i? gep = 2R to 1994 data from HI. 
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Figure 61: dependence of the dijet photoproduction cross section at Et = UGeV. We compare 

our NLO prediction with GRV and GS96 photon parton densities to 1994 data from HI using a fixed 
cone algorithm. 
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Figure 62: x^ dependence of the dijet photoproduction cross section at Et = 13GeV. We compare 
our NLO prediction with GRV and GS96 photon parton densities to 1994 data from HI using a fixed 
cone algorithm. 
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6 Numerical Results for Photon-Photon Scattering 



In this section we report on numerical results for inclusive one-jet and two-jet cross sections in 77 
collisions. We use the analytical results that have been calculated in leading order in section 3 and 
in next-to-leading order in section 4. The calculation proceeds as in the photoproduction case except 
that we now have three contributions: (a) the direct contribution, where both virtual photons interact 
directly with the quarks, (b) the single-resolved contribution, where only one of the photons interacts 
directly and the other photon is resolved, and (c) the double-resolved contribution, where both photons 
are resolved. It is clear that (b) and (c) have their analogs in the photoproduction case. For case 
(a) the analytical results are given in section 3.4 (Born matrix element), in section 4.1.3 (virtual 
corrections), and in section 4.2.9 (real corrections). 

Three partons appear in the final state of all three contributions when we include the NLO correc- 
tions. Two of these partons are recombined if they obey the cluster or the jet-cone condition. For the 
77 process we shall use only the Snowmass jet algorithm already introduced in section 2.4 with R = 1 
and R sep = 2. This algorithm was applied in the analysis of the TOPAZ, AMY [p, and OPAL Q 
data. However, there is no problem to introduce other jet definitions in the calculations, as we have 
done it for the jp case. Before we compare our results with recent data from the OPAL collaboration 
[ |9~i| ] obtained at LEP2, we shall present in the next section some general results and discuss some tests 
of the numerical program. 



6.1 Tests and General Results 

Similar to the "fp case, we have checked that our results are independent of the y-cut when the analytic 
results are added to the numerical results of the 2 — > 3 processes. This has been studied for all 2 — > 3 



processes separately. As an example, we show in figure |63| the y-dependence of the complete double- 
resolved two-jet cross section for LEP1 kinematic conditions and at T)\= f]i = 0> Et = 10 GeV. As 
to be expected, the 2-particle contribution is negative and decreases with decreasing y whereas the 
3-particle contribution is positive and shows the opposite behavior. The sum of both contributions is 
independent in the range 10" 4 < y < 10~ 2 . For lar ger y-cuts the independence breaks down because 
we neglected terms of order y in the analytic contributions. The slight decrease for y < 10~ 4 is caused 
by insufficient accuracy in the numerical integrations. This could be improved with more CPU time. 



For the results in figures p3 to 65, we use the LEP1 center-of-mass energy ^/s = 90 GeV. The 



photon spectrum is calculated from the formula 



n 



F 7/e (x) = - 



l + (l-x) 2 E 2 6 2 C {1 - x) 2 + m? e x 2 ( m 2 x 1 

In ~z ^ -\- 2(1 x) 



m 2 x 2 \E 2 9 2 (l — x) 2 + m 2 x 2 x) 

(6.1) 

where E is the beam energy and 9 C = 3°. The parton distributions are computed with the NLO 
parametrization of GRV [33] transformed to the MS subtraction scheme. For all scales we set /j, = 



M a = M b = E T and put N f = 5. AjyL = 130 MeV corresponding to A^L = 200 MeV as fixed in 

the GRV distributions. Inclusive single-jet cross sections have been calculated earlier J/], |95| . In these 
calculations, the double-resolved contribution has been obtained with a different method for canceling 



infrared and collinear singularities [73] which, however, can be applied only to the computation of this 
particular cross section. These results were compared to TOPAZ and AMY ||] jet production data 
and good agreement was found pjj. 

The fact that a different method for calculation of the double-resolved contribution to the inclusive 
one-jet cross section was available was utilized to test the new results based on the y-cut slicing method. 
Good agreement between the two independent methods was achieved for all double-resolved processes 
(<?<?', qq\ Q<1, qq, qg, 99 -> one jet) separately. 

In the following figures we show results for d 3 a"/d-E/rdr7id772 as a function of Et for special 7] values 
and as function of 772 for Et = 5 GeV and r\\ = using LEP1 kinematics. In figure ^3 we show the 
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Figure 63: Dependence of the inclusive dijet cross section on the y-cut, the boundary between analytical 
and numerical integration. Only the double-resolved part is shown. 



Et distribution for ryi = r/2 = in LO and NLO. We plotted the direct, the single-resolved, and the 
double-resolved contributions separately and show the sum of all three contributions. From the Et 
spectrum it is already visible that the direct contribution is reduced by the NLO corrections whereas 
for the resolved contribution the NLO corrections are positive. The correction in the double-resolved 
contribution is near 100%. This is seen more clearly in the lower part of figure 64, where we plotted 
the fc-factor as a function of Et for the three contributions and the sum. We remark that the LO 
curves are calculated with the two-loop a s and the same NLO parton distributions of the photon as 
in the NLO calculation. For a genuine LO prediction, one would choose the one-loop a s formula and 
LO parton distributions. This would change the fc-factors. In figure 34, the fc-factors show just the 
influence of the NLO corrections to the parton-parton scattering cross sections. The upper part of 
figure |64] also shows how the contributions of the three components (direct, single-resolved, and double- 
resolved) sum up to the full inclusive two-jet cross section. For relatively small Et, this cross section 
is dominated by the two resolved components. In this region we have a positive NLO correction and 
a fc-factor larger than 1. In the medium Et range and for large Et, the direct component dominates, 
the net NLO correction is negative, and the fe-factor is less than 1. It is clear that the importance of 
the resolved parts for larger Et increases with increasing center-of-mass energy. 

In figure 65, we show the rj2 distribution for 771 = and Et = 5 GeV, again for the three components 
separately and their sum in LO and in NLO. The ^-factors corresponding to the upper part of figure 
|65| are shown in the lower part. The three components show a very distinct behavior as a function of 
r/2. The single-resolved contribution does not have a plateau as broad as the other two components. 
So, by measuring in different 7/2 regions, it might be possible to enhance one or two of the three 
contributions. 

An important test of our calculations is the compensation of the factorization scale dependence 
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Figure 64: Direct, single-resolved, and double-resolved contributions and their sum in LO and NLO for 
the inclusive two- jet cross section at LEP1. Upper figure: Et -spectrum, lower figure: correspondent 
k-factors. 
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Figure 65: Direct, single-resolved, and double-resolved contributions and their sum in LO and NLO 
for the inclusive two-jet cross section at LEP1. Upper figure: n-spectrum, lower figure: correspondent 
k-factors. 
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between the direct and the single-resolved components and between the single-resolved and the double- 
resolved components. That this compensation works is shown in figures ^ a) and b) for the inclusive 
two-jet cross section at Et = 10 GeV and rji = 772 = 0- F° r this test we applied LEP2 kinematics, 
i.e. y/s = 175 GeV and 9 C = 1.72° in the photon spectra. In figure |6^ a), the NLO direct and 
the LO single-resolved cross section, where only the upper photon 7 a is resolved, are plotted as a 
function of M a /ET- The dependence of the NLO direct contribution on M a is clearly visible. This is 
opposite to the dependence of the single-resolved contribution originating from the scale dependence 
of the parton distributions. As one can see, the sum is constant as a function of M a /ET- The same 
compensation occurs between the NLO single-resolved and the LO double-resolved cross section. In 
the single-resolved cross section, only the lower photon 7^ is resolved. Also in this case (figure ^ b) the 
sum of both contributions is rather independent of M a . Another important topic for the NLO theory 
is the question of the overall scale dependence and whether this is reduced in NLO as compared 
to the LO cross section. First we show the dependence on the renormalization scale \x alone with 
the factorization scales M a = = Et fixed. In figure |67] a) it is clearly visible that in NLO the 
dependence on \x is reduced as compared to LO when [i is varied between /i = Et/2 and 2Et- If all 
scales [i = M a = = M are equal and this common scale M varies in the same range for the LO 
and NLO cross section we obtain the curves in figure |67] b). We see that the dependence in LO is such 
that the cross section as a function of M increases which is opposite to the behavior as a function of 
the renormalization scale /i in figure |67| a) . This comes from the dependence on the factorization scale 
which is not compensated (see figures |66| a) and b)). In NLO we have the effect that the dependence 
on the factorization scale is reduced due to the presence of the NLO corrections, so that in NLO the 
cross section as a function of M is fairly constant, i.e. it decreases only slightly with increasing M/Et- 
In conclusion we can say that the NLO cross section is nearly independent of the scales and therefore 
presents a much more solid prediction than the LO cross section. 

When high statistics data become available, one could test the variation of the jet cross section 
with changing cone radius R. This has been studied for the dijet cross section for Et = 10 GeV and 



f}i = V2 = 0. The result is displayed in figure 68 where the NLO cross section d 3 a/dETdrjidrj2 is shown 
as a function of R for R > 0.5. It increases as a function of R almost like a + blnR and somewhat 
stronger for R > 1. With our definition of the LO cross section the NLO result at R = 0.9 is equal to 
the LO cross section, so that the NLO corrections stay moderate if R varies between 0.7 and 1. 

6.2 Comparison with Experimental Data 

The first measurements of jet production in 77 reactions were done by the TOPAZ and AMY |l|] 
collaborations at TRISTAN. They presented data for the inclusive one- and two-jet cross sections as 
a function of Et in the range 2.5 GeV < Et < 8 GeV with different ranges of rapidity in the two 
experiments. These data were compared to predictions based on the theoretical work presented here 
in |J. Good agreement between the data of both experiments and the theoretical results was achieved 
for the one- and two-jet cross sections. The distribution in Et and also the absolute normalization 
agreed when the parton distributions were described by the GRV set. At the energy of the TRISTAN 
collider (y/s = 58 GeV), the main contribution to the jet cross sections comes from the direct process. 
However, the data could not be reproduced by the direct component alone. The resolved contributions 
were necessary to obtain agreement over the whole Et range. 

In the fall of 1995 the LEP ring was operated at the center-of-mass energy of 133 GeV. During the 
short run period, the OPAL collaboration collected data for jet cross sections in 77 collisions. They 
measured the inclusive one-jet cross section integrated over \rj\ < 1.0 and the inclusive two-jet cross 
section for 1 771 1 , 1 772 1 < 1 as a function of Et in the range 3.0 GeV < Et < 14 GeV J2|. These data 
were compared to our calculated NLO one- and two-jet cross sections in [g, and good agreement 
between the measurements and the predictions was found. The calculations were done with Nf = 5, 



A^jg- = 130 MeV and the NLO GRV parton distributions. In [11], also the dependence of the two-jet 
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Figure 66: Dependence of the inclusive two-jet cross section on the factorization scale M a in the photon 
7 a . Upper figure: NLO direct and LO single-resolved photon (only ^ a resolved), lower figure: NLO 
single-resolved (only 7^ resolved) and LO double-resolved photon. 
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Figure 67: Dependence of the inclusive two-jet cross section on the renormalization scale (i (a) and 
on the variation of all scales (i = M a = M& = M (b). 
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Figure 68: Dependence of the inclusive two- jet cross section on the size of the jet cone R for LO and 
NLO. 



prediction on other parton distribution sets was studied with the result that the other NLO sets, GS 



[]37j ] and ACFGP |34j], lead to almost the same results and the OPAL data could not be used to rule out 
one of these sets. In these calculations, the photon spectra were described by eq. ( |6.lD with 9 C = 1.43° 
as in the experimental setup. The cone radius was R = 1 and R sep = 2. 

In the meantime the LEP energy was raised and the OPAL collaboration extended their measure- 
ments to the center-of-mass energies of yfs = 161 GeV and 172 GeV Q[. They presented one-jet 
cross sections for |r/| < 1 and two-jet cross sections for |77i|,|t/2| < 1 (\/s = 161 GeV) and two-jet 
cross sections for \r]i\, 1 772 1 < 2 (^/s = 161 GeV and 172 GeV). The latter cross section was compared 
to our predictions already in [94]. As an example, we show the experimental results of the OPAL 
collaboration from the yfs = 161 GeV run. In figures 6£ and ff(] the inclusive one- and two-jet cross 
section data as a function of Et for jets with | ^7 1 < 1 are compared to the NLO calculations with 
R = 1, i?scp = 2, and the NLO GRV parametrization for the photon. The parameters were 9 C = 33 

(5) 

mrad, A^ = 130 MeV, and y. = M a = M b = E T . The direct, single-resolved, double-resolved cross 
IVlo 

sections and their sum (full line) are shown separately. As we can see, the agreement between data 
and the predictions is good. The resolved cross sections dominate in the region Et < 7 GeV, whereas 
at high Et the direct cross section dominates. We emphasize that the inclusive two-jet cross section 
is measured using events with at least two jets. If an event contains more than two jets, only the 
two jets with the highest Et values are used. This definition of the inclusive two-jet cross section 
coincides with what is done in the theoretical calculation. The good agreement between measured and 
calculated cross section is remarkable, since the NLO calculation gives the cross sections for massless 
partons, which are combined to jets, whereas the experimental jet cross sections are measured for jets 
built of hadrons. The good agreement then tells us that in 77 jet production, parton-hadron duality 
is realized to a high degree and no major disturbance due to underlying event energy is present. This 
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Figure 69: The inclusive one-jet cross section as a function of Et with \rj\ < 1 compared to our NLO 
calculation. The direct, single-resolved, and double-resolved cross sections and the sum (full line) are 
shown separately. 

is in contrast to what we observed for ^p jet production at lower Et where there was disagreement at 
the larger rapidity values. 

As a last point we confront the recently published | cos#*| distribution of the OPAL collaboration 
with the theoretical predictions. These data were taken for two x® BS intervals to separate direct and 
resolved dominated contributions similar as it was done for jet production in ep scattering (see section 
5). The selection cuts were > 0.8 (direct-dominated) and x^ < 0.8 (resolved-dominated), where 

T 2 i Et e~ r>i 

< = Affair- < 6 - 2 ' 

Y 2 i Et e m 

* = %^ (6 ' 3) 

are in LO equal to the fractions of the photon energies of the upper and lower vertex entering the 
hard parton-parton scattering process. Thus, in LO the direct process has x+ = x~ = 1, whereas the 
double-resolved process occurs only for x+,x~ < 1. The direct cross section data for | cos#*| between 
and 0.85 for the two event classes x^ > 0.8 and x^ < 0.8 are exhibited in figure [71]. The data points 



are normalized to have an average value of 1 in the first three bins and are plotted at the center of 
the bins. The data with x^ > 0.8 show a small rise with increasing |cos#*|, whereas the data for 
x^ < 0.8 show a much stronger rise in |cos#*|, similar to the findings in jp scattering, where we 
pointed out how this qualitative behavior is related to the exchanges of the LO cross sections. The 
theoretical curves for the cross section in the two x^ ranges are also shown using the two sets for the 
parton distribution functions, NLO GRV and NLO GS96 []3§|] . In these calculations, |r/i|,|r/2| < 1, 
|^| < 1, where fj = ^(r/i + 7/2), and the invariant two-jet mass Mjj > 12 GeV, as in the experimental 



selection [94]. Mjj is calculated from the two jets with highest Et- All other parameters are as stated 
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Figure 70: The inclusive two-jet cross section as a function of Et with |??i|,|??2| < 1 compared to our 
NLO calculation. The direct, single-resolved, and double-resolved cross sections and the sum (full line) 
are shown separately. 



previously. The theoretical curves in figure 71 are normalized at |cos#*| = 0. We observe that the 
shape of the | cos 6*\ distributions in the two bins is nicely reproduced in comparison with the data. 
This means that the qualitative behavior known from LO arguments is still present in NLO and agrees 
with the experimental data. On the quantitative side, the theoretical cross sections seem to increase 
stronger towards |cos#*| = 1 than the data, both the direct and the double-resolved event sample, 
indicate. However, a new OPAL analysis shows that the quantitative agreement there is much better 
for the resolved sample, if the /cr-cluster algorithm is used instead of the cone algorithm. We remark 
that the data analysis is still going on, i.e. the experimental data are still preliminary and we have to 
wait for the final analysis before further conclusions can be drawn. 



7 Summary and Outlook 

In this work we have presented a complete next-to-leading order calculation of direct and resolved 
photoproduction of one and two jets in jp and 77 collisions. Photon-proton and photon-photon 
scattering were considered simultaneously since the two processes are intimately related to each other. 
The results are of great importance not only as a test of quantum chromo dynamics (QCD), but also 
for the measurement of the proton and photon parton densities currently performed at HERA and 
LEP. 

First, we embedded the perturbatively calculable photon-parton scattering in the experimentally 
observable electron-proton scattering process, using the Weizsacker- Williams approximation and uni- 
versal structure functions for the proton and the photon. For these structure functions, we chose recent 
next-to-leading order parametrizations from the CTEQ and GRV collaborations. Special emphasis was 
given to experimental and theoretical ambiguities in the Snowmass jet definition. 
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Figure 71: | cos#*| dependence of the dijet cross section for Mjj > YlGeV, \rji\, \rj2\ < 1 and \r}\ < 1. 
We compare our NLO prediction with GRV and GS96 photon parton densities for direct- enhanced and 
double-resolved enhanced regions to preliminary data from OPAL. 



The hard photon-parton and photon-photon scattering cross section was calculated in leading and 
in next-to-leading order. The analytical calculation included the tree-level Born matrix elements, 
the virtual corrections with one internal loop, and the real corrections with the radiation of a third 
particle in the initial or final states. We integrated the latter over singular regions of phase space 
up to an invariant mass cut y. This was done in d dimensions in order to regularize the soft and 
collinear divergencies. All infrared singularities canceled, as they must according to the Kinoshita- 
Lee-Nauenberg theorem. The ultraviolet poles in the virtual corrections and the collinear poles in the 
initial state corrections were absorbed into the Lagrangian and the structure functions, respectively. 

The cross sections proved to be independent of the technical y-cut and less dependent on the 
renormalization and factorization scales, when added to the regular three-body contributions and 
integrated numerically. Excellent agreement was found in inclusive single-jet predictions with the 
existing programs of Bodeker and Salesch. We extensively studied the direct, resolved, and complete 
(in the case of jp) and the direct, single-resolved, double-resolved, and complete (in the case of 77) 
one- and two-jet distributions in the transverse energies and rapidities of the observed jets as well as 
the dependence on the jet cone size R. Finally, we compared similar distributions to inclusive and 
dijet data from HI and ZEUS and OPAL and found good agreement in all cases. 

It is possible to extend the work presented here in a number of different ways. First, the formalism 
presented can easily be extended to the calculation of "fp and 77 collisions, where one of the photons 
has a virtuality larger than zero, although still small compared to Et [10C]. Second, other observables 
than those considered in this work will allow to test the theory further and make it easier to isolate 
the resolved contributions to obtain information on the parton distributions in the photon for various 
scales. Third, our program can be used to predict single and dijet cross sections in proton-proton and 
proton-antiproton scattering at the TEVATRON at Fermilab in Chicago or at the designated LHC at 
CERN in Geneva. 
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A Phase Space Integrals for Final State Singularities 



This appendix contains the formulas needed to integrate the real 2^3 matrix elements in sections 
4.2.2 and 4.2.3 over phase space regions with final state singularities. To this aim, we have defined in 
section 4.2.1 a measure 







which contains all variables associated with the phase space of the unobserved subsystem p 1 = p\ +p%. 
The integral over the azimuthal angle <p is trivial, as the matrix elements are independent of <fi. 
We give the results for the four generic types of integrals in the following: 

= / d^F^ , a t = I d/iF^y / i T\ TR 

J z' z' + ab J z' z' + a(l — b) 

= ^ + \ (-i - \ m a ) + ma - \ + \ a - u 2 (-^) - ^ + 0(e)> (A,) 

/ 2 (a) = Jd^(l-b)(l-e) 

= -j- £ + \ + \lny F + 0(e), (A.3) 

/3(a) = J d ^^7 

= -i + lny F + 0(e), (A.4) 



/ 4 (a) = Jdfi F ^(l-b + b 2 ) 



= -i^-r 8 +0(e). (A.5) 
Terms of order e have been omitted since they vanish in the limit d —> 4. 



B Phase Space Integrals for Initial State Singularities 

In this appendix, we calculate the formulae needed to integrate the real 2^3 matrix elements in 
sections 4.2.5 through 4.2.8 over phase space regions with initial state singularities. The integration 
measure was defined in section 4.2.4 and is given by 

x a 

The integration variables are associated with the phase space of the unobserved particle p^. The 
integral over the azimuthal angle <f> is again trivial, as the matrix elements do not depend on <p. The 
longitudinal momentum fraction z a is integrated over numerically, since the cross sections still have 
to be folded with the parton densities. These are contained in the functions g(z a ) below. 
We give the results for the four generic types of integrals in the following: 



(1 " Z a ) 
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V 



1 - 2<j 



1 1 , 1 i 2 

2?"27 lna+ 4 ln a + Y 



+ 5(1) 
22(a) = y dmgiza)^ 

23(a) = J dii ig (z a )^(l- z a )(l- e) 



— + ln j// 

e V z n 



+ 0{e), 



X a 



dz a 

— g(z a ) 

z a 



--(1 - Za) + (1 " ^a) + (1 " ^a) In ( Vl — 



+ 0(e), 



, , 1 / Z^ _ 2 +2 _ £Z 2\ X 



dZr, 



9{Za) 



Xa 



l z 2 a -2z a + 2 + z^ - 2z a + 2 ^ ( \-z a \ _ 2 l-zg 

£ z a z a \ z a / z a 



(B.2) 



(B.3) 



(B.4) 



+ 0(e).(B.5) 



Terms of order e have again been omitted. 
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